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ABSTRACT 

We report the detection of Deuterium absorption at redshift 2.525659 towards Q1243+3047. We 
describe improved methods to estimate the Deuterium to Hydrogen abundance ratio (D/H) in absorp- 
tion systems, including improved modelling of the continuum level, the Lya forest and the velocity 
structure of the absorption. Together with improved relative flux calibration, these methods give D/H 
= 2.42±o;2!xl0 -5 cm" 2 from our Kcck-I HIRES spectra of Q1243+3047, where the error is from the 
uncertainty in the shape of the continuum level and the amount of D absorption in a minor second 
component. The measured D/H is likely the primordial value because the [O/H] = — 2.79± 0.05. This 
absorption system has a neutral Hydrogen column density log Nhi = 19.73 ± 0.04 cm~ 2 , it shows five 
D lines and is mostly ionized. The best estimate of the primordial D/H is 2.78^°' 38 x 10~ 5 , from the log 
D/H values towards five QSOs. The dispersion in the five values is larger than we expect from their 
individual measurement errors and we suspect this is because some of these errors were underestimated. 
We observe a trend in D/H with log Nhi that we also suspect is spurious. The best value for D/H is 0. Go- 
smaller than we quoted in O'Meara et al. (2001) from three QSOs, and although we have more values, 
the error is similar because the dispersion is larger. In standard big bang nucleosynthesis (SBBN), the 
best D/H corresponds to a baryon-to-photon ratio r\ = 5.9 ± 0.5 x 10 -10 and gives precise predictions 
for the primordial abundances of the other light nuclei. We predict more He than is reported in most 
measurements, although not more than allowed by some estimates of the systematic errors. We predict 
a 3 He abundance very similar to that reported by Bania et al. (2002), and we predict 3-4 times more 
7 Li than is seen in halo stars. It is unclear if those stars could have destroyed this much of their 7 Li. 
The r\ value from D/H corresponds to a cosmological baryon density Q^h 2 = 0.0214± 0.0020 (9.3%) that 
agrees with values from the anisotropy of the Cosmic Microwave Background: f^/i 2 = 0.021 ±0.003 from 
the Netterfield et al. (2002) analysis of BOOMERANG data and Q, b h 2 = 0-022±g"gg| from the Pryke et 
al. (2002) analysis of the DASI results. 



Subject headings: quasars: absorption lines 
cosmology: observations 

1. INTRODUCTION 



quasars: individual (cso 0167 = Q1243+3047) 



It is well established that the light nuclei hydrogen (H) , 
deuterium (D), 3 He, 4 He and 7 Li are all made during big 
bang nucleosynthesis. The relative primordial abundances 
created in the Standard theory of big bang nucleosynthe- 
sis (SBBN) for these five nuclei depend on one parameter, 
the cosmological baryon-to-photon ratio, 77= n&/n 7 (Kolb 
& Turner 1990; Walker et al. 1991; Schramm & Turner 
1998; Nollett & Buries 2000; Olive, Steigman & Walker 
2000). A measurement of the ratio of any two primordial 
abundances gives 77. The three other primordial abun- 
dances are predicted once 77 is known and measurements 
of them test the theory. 

SBBN has now been validated in two main ways. First, 
it successfully accounts for measurements of the approx- 
imate relative primordial abundances of all these light 
nuclei (Boesgaard & Steigman 1985; Walker et al. 1991; 
Copi, Schramm & Turner 1995; Schramm 1998; Schramm 
& Turner 1998; Tytlcr et al. 2000 and references there 
in). Second, the baryon density required is roughly con- 



sistent with that measured in other ways, including the 
Lya forest at high redshifts (Rauch et al. 1997; Weinberg 
et al. 1997; Zhang et al. 1998; Hui et al. 2002), and the 
baryon fraction in clusters of galaxies (Babul & Katz 1993; 
Boute & Canizares 1996; Bludman 1998; Rines et al. 1999; 
Arnaud & Evrard 1999; Wu & Xue 2000; Sadat & Blan- 
chard 2001). Recently, the cosmic microwave background 
(CMB) has given the same baryon density to within 10%. 

We use the D I/H I ratio in QSO absorption line sys- 
tems to estimate 77. This method is attractive because the 
deuterium abundance is more sensitive to 77 than other 
light nuclei, the D/H ratio in QSO absorption systems is 
apparently the primordial value, and we do not apply any 
corrections for unobserved ions. 

We measure D/H by dividing the column densities that 
we measure for absorbing D I and H I atoms. Since D/H 
is low, we could see D in only absorption systems that 
have high H I column densities, at most one per QSO. We 
would like to measure the column densities to better than 
10% accuracy and ideally < 1% to better test early uni- 
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verse physics. Echclle spectrographs give ample spectral 
resolution, and integration times approaching 1 day give 
the required signal-to- noise ratio (S/N). Once D has been 
found in a well calibrated QSO spectrum, there are four 
main factors that limit the measurement accuracy: the 
continuum level, the Lya forest, the velocity structure of 
the absorber, and contamination by absorption other than 
D. 

Each of the factors that complicate the measurement of 
D/H are hard to model and quantify. We need to know the 
level of the unabsorbed continuum to measure the amount 
of absorption. The Lya and other weak emission lines 
make it hard to estimate the precise continuum shape, 
and stochastic Lya forest absorption only makes matters 
worse. The Lya forest H I absorption is random in wave- 
length and opacity, and is ubiquitous over the redshift in- 
tervals in which we can measure D/H from ground based 
telescopes. All the H and D lines are in the Lya forest 
portion of QSO spectra, which introduces confusion and 
in the worse case can mimic D absorption. The hardest 
problem is to ensure that we explore all possible veloc- 
ity structures that might explain a given spectrum. D/H 
values and their errors are often dictated by the choice of 
velocity structure. 

It appears likely that the measurement errors have been 
underestimated in at least one of the published D /H val- 
ues, because these errors are hard to calculate, and the 
dispersion in D/H measurements is larger than expected. 
In O'Meara et al. (2001) we found that the dispersion be- 
tween three D/H values that we had measured towards 
different QSOs was larger than expected from the quoted 
errors. We concluded that this dispersion was probably 
not real, and instead a result of our underestimation of 
the errors. Since then, Pettini & Bowen (2001) have found 
D in a fourth QSO, Q2206-199 and quote a D/H value 
which is much less than the previous measurements, while 
D'Odorico et al. (2001) and Levshakov et al. (2002) have 
presented two very different measurements of D/H towards 
a 5th QSO, Q0347-3819. These new measurements fur- 
ther increase the dispersion in the reported D/H values, 
and again we suspect that this is because the errors have 
been underestimated. 

In this paper we announce the detection of D in 
an absorption system at z = 2.526 in the spectrum of 
Q1243+3047. We describe how we have improved our ex- 
ploration of some factors which determine the measure- 
ment errors. In particular we describe improved calibra- 
tion of relative flux levels, our estimates of the continuum 
level, our modelling of the velocity structure, and our mod- 
elling of the Lya forest and blended lines. 

We present material in the following order: observations 
and reductions (§2), overview and velocity structure of the 
absorption system that shows D (§3), measurement of the 
D column density (§4), evidence that we are measuring D 
(§5) measurement of the H column density (§6), the metal 
abundance and ionization in the absorbing gas(§7), the 
D/H from Q1243+3047 (§8), and from all QSOs (§9), and 
the deduced cosmological parameters (§10). The summary 
ends with a discussion of the improvements in methods 
that we use here, and a list of the issues that remain. In 
the Appendices we discuss how we model the continuum, 
our optimizing code, error estimates from the covariance 



matrix and D/H values from other QSOs. 

2. OBSERVATIONS AND DATA REDUCTION 

We report the detection of D absorption in the QSO 
Case Stellar Object CSO 0167, which we call by its B1950 
coordinate name: Q1243+3047. This QSO was reported 
as a V = 17 blue object at B1950 RA 12h 43m 44.9s DEC 
+30d 47.9', with a 15" error, (equivalent to J2000 12h 46m 
10.9s +30d 31m 31.2s) in an objective prism photograph 
by Sanduleak & Pesch (1984), and identified as a QSO at 
emission line redshift z em = 2.56 by Everett & Wagoner 
(1995). This QSO is not well known. The NASA extra- 
galactic database (NED) gives no other primary references 
except for a 2MASS QSO search. 

We have spectra with three different spectral resolu- 
tions, from three spectrographs, summarized in Tables 1 
and 2, and shown in Figures 1, 2 and 3. 

We obtained five spectra of Q1243+3047 from the Kast 
double spectrograph on the Lick Observatory 3.1m Shane 
telescope. These low resolution spectra are used to flux 
calibrate HIRES spectra. We compare these spectra to 
check the flux calibration, and we can sum them to im- 
prove the S/N. All integrations were obtained using the 
d46 dichroic that splits the spectrum near 4600 A, the 830 
line/mm grism blazed at 3460 A for the blue side, and the 
1200 line/mm grating blazed at 5000 A for the red side. 

We have one intermediate resolution spectrum from the 
ESI spectrograph (Epps & Miller 1998; Bigelow & Nelson 
1998; Sheinis et al. 2000), mounted at the Cassegrain focus 
of Keck-II. This echellette covers from 3900 - 11,000 A in 
a single setting in ten overlapping orders. We also use this 
spectrum to flux calibrate the HIRES spectra, and to look 
for metal line absorption at wavelengths larger than those 
covered in our HIRES spectra. 

We have eight high resolution spectra from the HIRES 
spectrograph (Vogt et al. 1994) on the Keck-I 10m tele- 
scope. These were obtained using the C5 dekker that has a 
1.14" slit width and gives a FWHM of 8.0 km s" 1 sampled 
with 2.1 km s" 1 per pixel on the original HIRES Tektronix 
2048x2048 CCD. To improve the relative flux levels along 
each spectrum, all observations were taken with the spec- 
trograph slits aligned with the local vertical. All but one of 
these spectra cover down to wavelengths below the Lyman 
limit break in the absorption system that shows D, while 
the exception extends to larger wavelengths to cover some 
metal lines. The S/N per pixel in the summed HIRES 
spectra (Table 2) increases linearly with wavelength up to 
4200 A, and then rises faster in the Lya emission line, 
reaching 105 on cither side of the Lya absorber of the 
D/H system. The S/N in the center of an echelle order is 
approximately 1.4 times greater than that at the ends. 

We describe the data reduction and flux calibration in 
detail in Suzuki et al. (2003). The error in the wavelength 
solution for the HIRES spectra is at least 0.05 - 0.1 pixels, 
or 0.1 - 0.2 km s _1 , and may be as large as 1 - 2 km s _1 . 
Our analysis of other similar spectra (Levshakov, Tytler & 
Buries 2000) revealed errors of order 1 km s^ 1 . We placed 
all the HIRES integrations on the same logarithmic wave- 
length scale, which has a constant velocity of 2.1 km s _1 
per pixel, similar to the original pixel size, and we shifted 
the Kast and ESI spectra to put them on a similar wave- 
length scale. All spectra were also converted to vacuum 



3 



wavelengths and shifted to the solar rest frame. We ap- 
plied a relative flux calibration to the HIRES spectra using 
the Kast and ESI spectra to transfer the flux information 
from standard stars to the QSO spectra. 

3. THE Z = 2.526 SYSTEM THAT SHOWS D 

In this section, we discuss the spectrum of Qf 243+3047, 
we introduce the absorption system that shows the D, and 
we discuss its velocity structure. 

3.1. Overview of the Q 1243 +3 047 spectrum 

We observe D in the conspicuous absorption system at 
z = 2.526 that produces both a clear Lyman limit break 
and has the strongest Lya absorption line in the spectra, 
near 4285 A. We sec no flux below the limit that is near 
3210 A in the Kast and HIRES spectra. No flux is ex- 
pected because the H I column density log Nhi is clearly 
near ~ 19 cm -2 from the shape of the Lya line. The pres- 
ence of a strong O I line at this redshift also implies that 
log Nhi> 19 cm~ 2 . In addition to five D absorption lines, 
this system shows many Lyman series H I transitions, but 
only a few metal lines, because the metal abundances are 
approximately 0.001 solar in the main component (see §7). 
In comparison, HS 0105+1619 has half as much H I, but 
shows 21 metal lines, as its metal abundance is approxi- 
mately 14 times larger (O'Meara et al. 2001). 

The system at z = 2.526 shares some of the properties 
of both the Lyman limit systems (LLSs) and the Damped 
Lya systems (DLAs). It is a LLS because it is optically 
thick in the Lyman continuum, and it has 25% of the 
minimum log Nhi to qualify as a DLA (Turnshek & Rao 
2002). It shows O I like DLAs, but unlike DLAs (Wolfe & 
Prochaska 2000), it is mostly ionized and the different low 
ionization lines have different velocity profiles. There is lit- 
tle published work on systems with log Nhi= 18—20 cm~ 2 , 
(e.g. O'Meara et al. 2001; Peroux et al. 2002) 

The spectrum shows metal lines from unrelated red- 
shifts, none of which effect our D/H measurement. A 
weaker Lya absorption line near 3710 k(z ~ 2.052) ap- 
pears to have weak damping wings that imply a high Nhi ■ 

3.2. How we Measure D/H 

The D/H value is the ratio of two column densities, Ndi 
for the neutral Deuterium (D I) and Nhi- We must esti- 
mate the errors on each of these quantities and show that 
both measurements apply to the same gas. 

We measure the Nhi to high accuracy from the Lya 
line alone, because the shape of this line is dominated by 
damping wings that are sensitive to Nhi, and less sensi- 
tive to the velocity information. However, the Nhi value is 
more accurate and reliable when we use accurate velocity 
information, and hence we measure the Nhi last, after we 
have examined the velocity information. Unlike the Nhi, 
the Ndi that we measure depends strongly on the velocity 
information. 

We obtain the velocity information from the D I lines 
themselves, taking guidance from the high order H Lyman 
lines and the metal lines, especially O I. 

3.3. Velocity Structure of the z = 2.526 absorption 
system 



We now discuss the clues that metal line absorption 
give for the velocity structure. The system at z — 2.526 
shows many components, with different structure in differ- 
ent ions, but the bulk of the H I is in just two components 
that we see in absorption by O I. 

We expect that the H I and the D I should have identi- 
cal velocity distributions, except for the effects of thermal 
broadening. We expect the D/H ratio to be constant from 
component to component, because the metal abundances 
are too low for significant destruction of D. We also ex- 
pect that Ndi /Nhi = D/H because the D and H should 
have nearly identical ionization (Savin 2002). 

We would like to measure the velocity structure of the H 
and D separately, allowing for a direct comparison. How- 
ever, the H lines are all saturated, broad and blended, and 
give little information. Instead, we use the unsaturated 
metal lines, especially the O I lines, as guides. 

We do not know how closely the metal lines will trace 
the velocity structure of the H I and D I. We expect to 
find some H I wherever we find metals, but the metal lines 
can have different velocity distributions in detail because 
the ionization and metal abundance can vary from com- 
ponent to component, and perhaps with velocity inside a 
component (§5.1). 

In Figure 3, we present the regions of the spectrum 
where we expect metal line absorption. We observe strong 
absorption in only a few metal ions: O I, C II, and Si II. 
C III and Si III also show absorption, but are poorly con- 
strained since their lines appear saturated and may be 
highly contaminated by Lya forest absorption. We see 
weak C IV and Si IV absorption that has very different 
velocity structure from the low ionization metals. Al II 
is seen in the ESI spectrum, but is not covered by our 
HIRES data. Fe II 1608 is not detected at our S/N, and 
although some absorption appears near the expected posi- 
tions of N I, N II, N III, and Fe III, it is very weak, and can 
readily be accounted for as Lya forest. We also measure 
\ogNc ii* < 12.3 cm~ 2 in all components. 

The O I line suggests that two components will be 
needed to model the velocity distribution of the gas that 
shows the D. O I provides the best indication of the ve- 
locity distribution of the H I and D I absorption, because 
O I/H I is similar to O/H in gas of low ionization (O'Meara 
et al. 2001). The O I 1302 transition is in a high S/N re- 
gion of the spectrum well separated from other lines. This 
line is asymmetric, with extra absorption at larger wave- 
lengths. We fit the O I with the two components that we 
list in Table 3. 

In Figures 4 & 5 we show fits to the O I, C II, Si II, 
C IV, and Si IV lines. Voigt profile fitting usually does 
not produce unique results, and this is especially true for 
C II and Si II, which can be synthesized with a variety 
of components at different velocities that are all heavily 
blended. The C II and Si II fits listed in Table 3 and 
shown in Figures 4 & 5. Irrespective of the exact model 
used, the C II and Si II require gas at v ~ —40, ~ 0, and 
~ 95 km s , blended with additional absorption between 
16 and 35 km s _1 . 

The C IV and Si IV lines have a different velocity struc- 
ture from the low ionization transitions O I, C II and 
Si II, which means they primarily come from different 
gas. Both C IV and Si IV show absorption at v ~40, 
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100, and 200 km s -1 . There is also a component centered 
near v ~ 5 km s _1 for C IV and v ~ —10 km s _1 for Si IV. 

In summary, the metals indicate that H I may be found 
in components near: v = —40, 0, 13, 20, 40, 95 and 
200 km s _1 . We are confident that the O I marks the 
velocity of most of the H I. The low ionization lines of 
C II and Si II have most of their column density at these 
velocities, and O I is not seen at other velocities. 

There is no sign of any metal components near v ~ 
—81.6 km s" 1 , the velocity that would place H I lines at 
the position expected for the D I lines. This means that 
there does not need to be any H I that could contaminate 
the D lines, although there could be if the metal abun- 
dances were very low. 

3.4. Velocity Structure of the Gas that Gives D/H 

While the system shows many components, the D/H 
measurement will depend on just three. We expect to see 
D I in the two components that we see in the O I, C II and 
Si II, while the third component is an unfortunate blend 
that adds uncertainty but does not show D. 

• Component 1, near v = km s™ 1 , is the strongest 
component in O I, C II, Si II and D I. We will re- 
fer to the D I absorption in this component as D-l, 
with column density Ndi(D-I) and velocity disper- 
sion &(D-1) and similarly for the other ions. The 
parameters of the O-l and D-l components are well 
determined and they contain approximately 90% of 
the column density. We use the redshift of the O-l 
component, z (0-1)= 2.5256916, to define a refer- 
ence frame with radial velocity v — km s _1 . In 
this frame, we measure the velocity of the O-l as 
v = 0.00 ± 0.14 km s _1 , where the error is from the 
fit alone and ignores the poorly known wavelength 
scale error. 

• Component 2, near v = 13 km s _1 , is seen as an 
asymmetric extension of the O I that we call com- 
ponent 0-2. There is also some C II, Si II near 
this velocity. We call other associated lines D-2 and 
H-2. Neither the velocity nor the b- value of the O- 
1 component are well known, and they are corre- 
lated. Smaller 6-values are needed as the line cen- 
ter moves to larger wavelengths. If we have some 
additional velocity or 6-value information, then the 
second component is well determined. For example, 
if both O I components have the same b- value, this 
value is 6.78 km s^ 1 , component 2 is at 13.3 km s _1 , 
and its column density at 13% of the total. 

• Component 3, near v = —40 km s _1 , is an even 
weaker component in C II and Si II. We will refer 
to the H I absorption from this component as H-3, 
which is at v(H - 3) = -44.3 ± 7.0 km s" 1 , with 
log N m (3) = 15.90 ± 0.03 cm" 2 from the higher or- 
der Lyman lines. This log Nhi value is low enough 
that we would not expect to see O I or D I and 
neither is seen. H-3 effects D/H because the short 
wavelength side of the Lyman lines from H-3 acci- 
dentally blend with the long wavelength side of the 
D-2 lines, and hence some of the absorption near - 
70 km s" 1 can be explained by either H-3 or D-2. In 



Figure 6, we see that D-2 is strongly blended with 
D-l and H-3. 

There is a large amount of H I absorption, log Nhi— 
19.7 cm~ 2 , near components 1 and 2, but we have lit- 
tle information on the velocity distribution of this H I. In 
Figs. 7 & 8 we show the Lyman series lines. The spectra 
require H I Lyman lines corresponding to the metal line 
components near v = —40, 0, 95 and 200 km s" 1 . The 
spectra are also consistent with, but do not need, com- 
ponents at 13, 20 and 40km s -1 . The parameters that 
we find for these components are in Table 3. If there is 
H I absorption at 40 km s _1 with b > 10 km s _1 , the 
blue edge of the higher order Lyman lines show that it has 
log Nhi < 16.4 cm" 2 . Figure 8 shows there is no flux in 
the range —40 to 40 km s _1 in lines up to Ly-14, and in 
higher order lines the S/N is rather low, and the Lyman 
lines start to overlap. We can measure the shape and width 
of the higher order lines, but unfortunately this tells us lit- 
tle about the gas at v = — 13.3 km s _1 . The absorption 
on the negative velocity side of these lines is determined 
by the H I at -40 km s _1 with log Nhi ~ 16 cm -2 , and 
the positive velocity side of these lines may be influenced 
by log Nhi ~ 16 cm -2 absorption at 40 km s _1 . 

A single component at v = km s _1 with a b ~ 
17.7 km s _1 gives a fair fit to all of the H lines without 
separate lines at ±40 km s _1 . This b- value is an upper 
limit on the b- value for the H I in component 1, since we 
must use a lower b- value when we include separate lines at 
±40 km s _1 . We will see below that we can also obtain 
an excellent fit to the H lines using components H-l and 
H-2. 

4. MEASUREMENT OF THE D COLUMN DENSITY 

In this section we explore a variety of models that can 
account for the spectrum. We estimate the log Ndi from 
the models that fit the spectrum with the lowest \ 2 values 
(§4.1) and in §4.2 we discuss the many factors that effect 
the accuracy and reliability of the log Ndi value and our 
estimate of the error. 

To make a reliable estimate of the D/H ratio in this 
system, we must explore the full range of models that can 
explain the H and D absorption in the spectrum. Different 
models employ different assumptions about issues such as 
the velocity structure of the absorber, and whether D or 
H is making parts of the absorption. 

A large part of the uncertainty over the D/H comes from 
the log Ndi in the D-2 component. We can not predict the 
amount of D in D-2, because we do not know the fraction 
of the H I in this component, and the O/H need not be the 
same in the two components. The spectrum requires some 
absorption near -70 km s _1 that can be either D-2 or H-3. 
We must now determine how much of this absorption is 
from D-2. 

Most of the information on the D is in transitions Ly-2 
to 8, which we show in Figure 7. In each transition we 
see absorption that has all the expected characteristics of 
component D-l associated with H-l and O-l. In Ly-4 to 8 
we see a resolved D line that is well separated from the H 
absorption, while in Ly-2 and 3 the D absorption is fully 
blended with the H lines, as we expected because the H line 
is wider. The velocity structure, central velocity, 6-value 
and column density are all well determined from Ly-4 to 
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8, especially since the last few lines are not saturated. The 
absorption is coming from D in gas with velocity near zero. 
The D-l lines are all very narrow, and they can all be fit 
by a single component. 

However, we have less information on the D-2 compo- 
nent, because it is fully blended with the D-l on its short 
wavelength side, and H-3 on the other side. 

In Figure 9 we show the b and log Ndi values that we 
obtain when we fit each of the D lines individually with a 
single component. While there is general agreement, the 
dispersion is clearly larger than expected for the errors 
that we show. These errors come from the covariance ma- 
trix of the x 2 optimization and we explain in the Appendix 
§14.1 why such errors are often too small. 

Two of the D-l lines are contaminated by other lines. 
The Ly-5 D line contains a second narrow Ly-3 line at 
z = 2.3987515 that has log N H i 14.459 cm" 2 and b = 
28.81 km s _1 . This contaminant is well determined by 
its other Ly lines and we fit it when we fit the Ly-5 D 
line. The D Ly-3 line includes a slight contribution from 
a Ly/3 line with z = 2.342177, b = 35.6 km s" 1 and 
log N H i = 13.81 cm" 2 . 

4.1. Grid Search to Find log Ndi 

In this section we estimate the log Ndi value and its er- 
ror. Wc have explored various models that might explain 
the absorption at all wavelengths relevant to deuterium. 
When we calculated the % 2 difference between the simu- 
lated spectra from the models and the data we found that 
there was a well defined set of parameters that give excel- 
lent x 2 values. We obtained log Ndi and its error from the 
set of models that had x 2 < Xmin + 1- 

We made a six dimensional grid of models, one model 
for each combination of the allowed parameter values. We 
specified these models using six parameters that we vary, 
because they might change the D/H, and other parame- 
ters that we hold constant. We made simulated spectra 
for 15,750,000 different models, one for each possible com- 
bination of the parameter values. We allowed each of the 
six parameters to have every discrete value given by the 
(minimum value, maximum value, and the number of uni- 
formly spaced steps), as follows: 

• log N D i (total) (15.07 cm- 2 , 15.2, 15). 

• Ndi (D-2) / N D i (total) (0.01, 0.35, 35). 

• v sep , the velocity separation of the two D compo- 
nents, (5 km s , 20, 20). The maximum value lets 
the D absorption cover the full range of the O I ab- 
sorption, while we found that log Ndi (D-2) was in- 
significant whenever the separation was < 5 km s . 

• 6(D-1) (9 km s-\ 9.48, 10). 

• 6(D-2) (4 km s _1 , 15, 10). The upper limit is guided 
by the largest value for 6(11-2) and we discuss the 
lower limit in §4.2.4. 

• 6(H-3) (15.7 km s"\ 19, 15). 

We found that the log Ndi (D-l) varies very little, as ex- 
pected, because its lines are strong and clear. The main 
variation in log Ndi (total) comes from log Ndi (D-2), and 



hence the first two parameters are strongly correlated. All 
except the first and last parameters in the list describe the 
velocity structure of the D, while the last parameters varies 
the amount of H absorption near the D-2 lines. In most 
cases the ranges and the step sizes were chosen by trial and 
error to sample all models that can give fits with low x 2 ■ 
Exploring the range of parameters in this grid search was 
a computationally intensive task that took over a month 
and required tens of runs, many of which took over 24 
hours on a 1.5Ghz personal computer. We used a similar 
grid search, with only three parameters, in Kirkman et al. 
2001. 

Several other parameters of the model were not varied, 
because we believed they would have less effect on D/H, 
and the computations take approximately ten times longer 
for each additional parameter. The fixed parameters in- 
clude: 

• The wavelength ranges included in the fit: -40 to 
-120 km s _1 for Ly-2 to Ly-8. 

• v (D-l) = —2.8 km s _1 , which is well constrained by 
the spectrum and is also used as the value for i>(H- 
1). We discuss why this differs from the O-l velocity 
in §5.1. 

• The specification of the main H I, including its 
log Nhi (total) and its two components, H-l and H- 
2. We use the values given in Tabic 3, some which 
come from the fits that we will describe in §6. We 
can fix the main H I because the high order Lyman 
lines of the H I are all strongly saturated and they 
do not absorb near the D lines. 

• log Nhi (H-3) is constrained because H-3 does not 
absorb in the higher order Ly- lines (Figs. 6 & 8) 
and its precise Nhi value does not effect the log Ndi- 

• u(H-3) = -44 km s" 1 . Wc found that 6(H-3) is ef- 
fectively degenerate with v(H-3), and hence cither 
parameter suffices to explore the range of accept- 
able models. We changed w(H-3) from -50 km s -1 
to -30 km s _1 and found that the &(H-3) was larger 
in compensation, but there was no change in the 
log Ndi (total) that gave acceptable x 2 values. We 
believe that this is a sufficient range to explore for 
v(H-3) because a metal component is centered near 
-41 km s -1 . 

• The continuum shape and height. In Figure 10 we 
show the continuum that we used near the Ly- limit, 
while in Figs. 6, 7, 8 & 11 wc show the fit to the 
spectrum using this continuum and absorption by 
obvious lines. We found that the bumps in the con- 
tinuum are required to fit the spectrum. We do not 
claim that the QSO continuum itself is this bumpy, 
rather that the combination of the QSO contin- 
uum, emission lines, broadly distributed and weak 
blended absorption, and flux calibration errors com- 
bine to make this the effective continuum. Had we 
chosen a smooth or flat continuum across Fig. 10 
then we would have also required extra broad ab- 
sorption distributed like the dips in the continuum 
in Fig. 10, or emission near 923 A in the rest frame 
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of the QSO, to obtain fits similar to those seen in 
Fig. 11. 

• The absorption in the vicinity of D that is H I at 
different redshifts. 

In §4.2 we show that the last two factors are related and 
have little effect on log Ndi- 

In Figure 12 we show contours of constant x 2 as a func- 
tion of the total Ndi against each of the other five param- 
eters that we adjust. These plots show projections of the 
minimum x 2 values onto the planes specified by the axes, 
and hence the values of the four other parameters may vary 
slightly as we move along a particular contour. We found 
these minimum x 2 values by stepping each adjustable pa- 
rameter over its full range, to explore all combinations of 
all parameters. 

Our best estimates for the parameters are the values 
that together give the minimum % 2 value, Xmin- These 
values arc all contained within the innermost contour on 
each panel. The ma confidence intervals for the parame- 
ters are the ranges allowed by the successive contours on 
the plots, that we have placed at values of Xmin + m2 - 
Hence the la ranges are the smallest and largest values 
of the parameters that are allowed by the innermost con- 
tour, at Xmin + 1- These critical x 2 values apply because 
we have varied all the parameters to find the minimum x 2 
value for each log Ndi (total). We have only one effective 
degree of freedom, the log Ndi (total). The distribution of 
the x 2 statistic as a function of log Ndi (total) is then ex- 
pected to be approximately the Xi distribution function, 
with one degree of freedom. We shall refer to this method 
of estimating the errors on a parameter as the Sx 2 method. 

In Figure 13 we show the minimum x 2 that we found 
for each log Ndi ■ The critical x 2 values are the same as in 
Figure 12, for the same reason. Our best estimate for the 
total D in components 1 and 2, the column density that 
we will use to get D/H, is 

log N BI (total) = 15.113±°;^ 2 cm" 2 , (1) 

which is a la range of +10% and —6%. This is the 
log Ndi (total) value that gives the lowest x 2 when com- 
bined with the following values for the other parame- 
ters, all of which arc collected in Table 3): 6(D-1) = 
9.2±0.2km s _1 , 6(D-2) = 4 km s~\ v sep = 12±| 5 km s" 1 , 
N D i(D-2) / Ndi (total) = 0.12+°;^, and hence log N D i(D- 
2) = 14.191 cm" 2 . 

The Xmin = 271.4 is reasonable for 274 pixels. We have 
used 3 degrees of freedom fitting the parameters of each of 
the two D lines and 2 parameters for the H-3 lines. The 
continuum, other contaminating lines, and Nhi(H-3) are 
well constrained by pixels outside the region used in the 
grid search, while v(H-3) is degenerate with 6(11-3). There 
are then approximately 274-8 degrees of freedom. The 
probability xlee > 271 - 3 = °- 4 0. 

We have a simple explanation for the range of 
Ndi (total). The x 2 distribution is asymmetric, rising more 
slowly as log Ndi (total) increases from the best value, be- 
cause of the contribution from D-2. The log Ndi(D-I) is 
well determined and varies little for all models that give a 
small x 2 while the log Ndi (D-2) varies from zero to a well 
determined maximum. Figure 12 shows an anti-correlation 
between log Ndi and 6(11-3), which arises because the ab- 



sorption near -70 km s _1 can be explained by various com- 
binations of D-2 and H-3. As log Ndi (D-2) decreases, the 
6(H-3) can increase to give models with a similar shape 
and x 2 . The lower bound on log Ndi comes from our as- 
sumption that the absorption near v = is D, which we 
discuss in section 5. The lower bound on log Ndi applies 
when log Ndi (D-2) is insignificant and all the absorption 
in D-l is D. The upper bound on log Ndi is met when D-2 
absorbs photons that we see at velocities ~ 70 km s -1 . 
We find a continuous range of acceptable models between 
these two bounds. 

The fit shown in Figures 8, 6, 11, & 7 is that with the 
minimum x 2 ■ Other fits with x 2 values larger by several 
look indistinguishable in such plots. 

We found that we could also find good fits to the D 
lines using the optimization methods that we describe in 
§6. Before we conducted the grid search, our best fit had 
a reasonable x 2 = Xmin + After we completed this 
grid search, we confirmed that the optimization methods 
can find the same best fit, but only when we limit the 
parameter to narrow ranges centered on the best fit. 

4.2. Discussion of the log iVoi Value and its Error 

We now discuss a few of the many factors that can ef- 
fect the log Ndi value and its error. In each case we have 
made some reasonable choice for the relevant parameters 
when we conducted the grid search. Here we discuss some 
of these other options that we did not quantitatively ex- 
plored in the grid search. In all cases we suspect that 
these effects are smaller than the error that we quote on 
the log Ndi, except for missing D components that might 
be comparable to the error. Some of these factors could 
appear as random or systematic errors, depending on de- 
tails such as potential bias. An error is random if its effect 
on log Ndi is such that the mean log Ndi from many spec- 
tra and model fits would converge on the true value. 

4.2.1. Flux Calibration 

Flux calibration is a possible cause of the bumpy contin- 
uum in Fig. 10. Four different HIRES orders contribute at 
A < 3400 A. The wavelength scale of the bumps are con- 
sistent with the 10 A correlation length that we expect for 
flux calibration errors in this part of the spectrum. The 
amplitude is also consistent, except for the largest bump 
near 3290 A that is 6 - 12% in amplitude, two - three times 
larger than the expected calibration errors. There are two 
components to the flux calibration error, the S/N of the 
HIRES integrations, and the S/N of the Kast spectra. We 
find that the total error in the calibration of the summed 
HIRES spectrum is approximately a 4% near 3290 A, and 
8% near 3250 A where the S/N is lower. 

The method we used to fit the continuum near the D 
line should have corrected for the errors in the flux cali- 
bration without adding to the error on the log Ndi- The 
continuum fitting will also have corrected for any other 
effects that are correlated over > 10 A. 

4.2.2. Continuum Level 

We find that the continuum level is not a significant 
part of the error in the log Ndi value. If we multiply the 
adopted continuum by 1.05, the result is clearly too high 
for all the D line, while if we multiply by 0.95 the result 
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is too low, except perhaps for Ly-3. This means that the 
error in the local continuum is < 5%, and since we fit the 
continuum independently for each of the 5 main D lines, 
the corresponding error on the log Ndi is < 5%/y/5, which 
we will ignore. 

4.2.3. Velocity Error 

We took v(D — 1) = —2.8 km s^ 1 and we did not vary 
this value, to account for the error of 0.6 km s" 1 from the 
fit to D-l or the wavelength scale error. If the D-l were not 
centered at this velocity, it would have a smaller log Ndi 
and some additional narrow D or H absorption would be 
needed to account for all the absorption. This does not 
seem a likely way to change the log Ndi value. 

4.2.4. Could b(D-2) be very Low? 

The grid search hints that 6(D-2) could be very small. 
The model with Xmin had b(D-2) = 4 km s _1 , but with a 
la range extending to 12 km s" 1 . The minimum b value 
that we considered in the grid search was 4 km s _1 be- 
cause the instrumental resolution is similar, and hence 
b < 4 km s _1 has little effect on the spectrum. 

Small values for 6(D-2) do not allow very large 
log Ndi(D-2). The spectrum places an upper limit on the 
opacity near -70 km s , and a weaker limit on the equiv- 
alent width. We can then trade Ndi(D-2) against 6(D-2) 
to find acceptable fits for a wide range of 6(D-2). Solutions 
with low 6(D-2) tend to have low Ndi(D-2), the opposite 
of the situation for lines with constant equivalent width. 
Fortunately the optical depth prohibits large log Ndi (D- 
2), and hence the uncertainty in the 6(D-2) value has little 
effect on log Ndi(D-2), as seen in the third panel of Figure 
12 . 

Low b(D-2) values require temperatures that are un- 
likely. If 6(D-2) < 4 km s^ 1 , the absorbing gas would have 
T < 1900 K, and it should be neutral if ionized by inter- 
galactic radiation. However, the C II/O I and Si II/O I 
ratios appear to be higher near to component 2 than in 
component 1 , which implies that component 2 is more ion- 
ized, and less likely to be neutral (§7.2). 

4.2.5. Velocity Structure 

When we calculated the maximum log Ndi in D-2, we 
assumed that both D-2 and H-3 were Voigt profiles. More 
complex velocity distributions could give different column 
densities, a topic that has been explored by Sergei Lev- 
shakov and collaborators (Levshakov et al. 2003). This is 
unlikely to be a major effect for H-3, which is wide enough 
(Fig. 12) that its profile is dominated by thermal motions, 
and hence should be insensitive to the detailed velocity 
distribution. The same argument wil apply to D-2 if its 
b- value is close to the la upper limit of 12 km s^ 1 , but not 
if 6(D-2) is near 4 km s _1 . 

4.2.6. Missing D Components 

The D/H would be systematically low if we have in- 
cluded velocity components in log Nhi (total) that are not 
detected in D I because the D I is hidden by H I absorp- 
tion. 

We see no evidence that we have missed D components, 
and we do not know how to estimate the chance that they 



exist, and hence we will not include them in the error 
on log Ndi- However, they could exist and have a major 
effect on log Ndi, even doubling the D/H if their metal 
abundances are very low. 

We need to examine the D and H lines in detail to es- 
tablish the velocities that could contain missing D com- 
ponents. Fig. 8 shows that the H I included in the total 
log Nhi ah lies well within the range -40 to 40 km s _1 . At 
other velocities, log Nhi< 17 cm -2 because see flux at the 
expected positions of one or more of the higher order Ly- 
man lines. Figs. 7 & 11 show that there are no significant 
D lines from components at -40 km s^ 1 < v < km s _1 , 
except for D-l. At positive velocities there are opportuni- 
ties to hide D. Component 2 was hard to detect, some of 
component H-3 could be D, and there might also be sig- 
nificant D between the two, at 13 < v < 40 km s^ 1 , where 
we see C II, Si II, C IV and Si IV. 

The H I lines alone allow that a major portion of the to- 
tal log Nhi is at velocities approximately 16 - 25 km s _1 . 
H I lines alone do not give the v of the main H I to within 
approximately 25 km s _1 because the core of the Lya con- 
tains many other components, and the higher order Lyman 
lines are blended on both sides. However, the C II, S II 
and especially the O I all indicate that most of the H I is 
in component H-l. The H I in components other than 1 
and 2 is probably < 10% of the total because we do not 
see O I at other velocities. If D/H is a constant, we expect 
that any missing D is also less than 10% of the total. We 
would have seen the O I from any component with more 
than 10% of the total log N H i provided it had [O/H] > -3. 

We do not know whether the absence of O I and other 
metals is sufficient to rule out missing components that 
contain > 10% of the total column density. Abundances 

< — 3 are very rare for absorption systems as a whole. Not 
one of 34 DLAs listed by Pettini et al. (1997) had [Zn/H] 

< —2, although several were upper limits and some sys- 
tems with lower abundances are known (e.g. Fan & Tytler 
1994). However, we know little about the abundances in 
components of systems. 

Although H-3 might contain as much optical depth in 
D as D-l, H-3 is predominantly H and not D because 
a distinct component of C II and Si II is centered near 
v = -40 km s" 1 (Table 3, Figure 4). 

4.2.7. Contamination of the D-l component by H 

In §5 we explain why it is unlikely that the D-l lines are 
contaminated by H, other than the blends that we men- 
tioned above with D Ly-3 and 5. Several of the D lines 
seem relatively free of close contamination and hence we 
do not expect that the way in which we fitted the contam- 
inants to D Ly-3 and 5 is a significant source of error. 

4.2.8. Validity of the Model and the x 2 Constraint 

The error we quote for log Ndi comes from the range of 
fits with x 2 < X 2 m in + 1- The log Ndi and the error will 
be unreliable if: (1) the model or the parameters with the 
lowest x 2 are unrealistic, or physically unacceptable; (2) 
we did not fully explore the parameter space of the model; 
(2) we did not explore a wide enough range of models, for 
example allowing D in enough components or at enough 
velocities. 

We believe that the log Ndi value and its error are rea- 
sonable. We explored a relatively large range of possible 
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models using coarse grids before we settled on the fine grid 
presented above. However, we may have misunderstood 
the velocity structure because the S/N is low near the Ly- 
man limit and we did not explore more complex models 
because they would have taken too long to compute. 

5. THE ABSORPTION NEAR THE D-l POSITION IS 
MOSTLY D 

Several lines of argument imply that the absorption that 
we identify as D-l is mostly D. The evidence is convincing, 
but not as strong as for some other QSOs. 

Ideally, we would compare the velocity structures of 
the D I and H I lines, which would be identical if we 
have correctly identified the transitions and there is min- 
imal contamination. This comparison is not possible for 
Q 1243+3047 because the velocity structure of the H I is 
not observed, and hence we are less certain that we have 
seen D than we were in other QSOs. 

The only reasonable identifications for D-l are D and 
H, because D-l shows a Lyman series. 

If the D-l absorption is H, 6(D-1) implies a very low 
temperature of < 5130 K. This temperature is much too 
low for low density gas photoionization by the intergalactic 
ionizing UV radiation. It is certainly too low for the Lya 
forest and perhaps for components with log Nhi= 15 cm~ 2 
in DLAs and LLS, unless they arc shielded from the radi- 
ation by H I. There are no metals at the redshift of H at 
the position of D-l, but this is not significant, since metals 
would not be seen if their abundances were low. 

There are three main reasons why D-l is D: its velocity, 
line width, and column density are all close to the values 
that we expect. 

The column density of the absorption near D-l could 
have been orders of magnitude different from the measured 
value. D-l is likely to be D because it has a log Ndi that 
gives D/H similar to the values found toward four other 
QSOs. The converse is not true: we would not reject a D 
line that gives a different D /H value if it otherwise appears 
to be D. 

5.1. Velocity Agreement 

The wavelengths of the centers of the D, O and H lines 
indicate that the D-l absorption is D, but they also suggest 
that O /H varies with velocity. 

In the frame of the main component of the 0-1, D-l is 
at —2.8 ± 0.6 km s _1 , where the error is from the line fit 
alone. This is close enough to show that D-l is mostly D, 
since the D-l and 0-1 lines overlap in velocity. 

The O-l, D-l velocity difference of 1.3 pixels is larger 
than wc expect from measurement errors, but we have 
been surprised in the past by wavelength scale errors in 
HIRES spectra. The error might be with the single O I 
line, rather than the D lines, since the several D lines ap- 
pear at the same velocity (Figs. 6 & 7). 

The H I velocity is consistent with that of the D I, but 
it is not known well enough to show that the D-l line is 
D. 

5.1.1. Does O/H Increase with Velocity? 

The O I might have a different velocity structure from 
the H I and the D I if O/H varies and is correlated with 
velocity. For example, we can imagine that all of the H I, 



D I and O I come from a single component. This compo- 
nent might have a Gaussian distribution of velocities. The 
0-1 component would be centered at +2.8 km s _1 in the 
frame of the D-l because O/H is larger at larger velocities, 
and the 0-2 component might arise from the O/H gradi- 
ent alone. In this case we should model the H and D with 
a single component. We do not favor this model, because 
C II and Si II show components 1 & 2. However, we can 
fit to the H I lines with a single component (§3.4), and 
this can have a 6-value that is consistent with the D and 
O lines (§5.2). If we have a single component the log Nhi 
is unchanged (§6), while the D/H is lower because log Ndi 
tends to be lower when we loose D-2 (§4.1). 

5.1.2. Metal Abundance Gradients in Other Absorbers 

There is almost no published information on the varia- 
tion of metal abundances or metal line kinematics in com- 
ponents of absorption systems with log Nhi— 19 cm~ 2 , 
especially on velocity scales of a few km s _1 . The closest 
analogue is HS 0105+1619 where we found that ions C II, 
N II, Si II and Fc II did not come from the same gas as 
the H I, D I and O I, because the line centers differed by 1 
- 2 km s _1 and the lines widths were different (O'Meara 
et al. 2001). Toward absorption systems with much lower 
Nhi we have found that the abundances can vary greatly 
between components with similar velocities and relatively 
high ionization. Toward PKS 1937-1009 Tytler, Fan & 
Buries (1998) found that two components separated by 
15 km s _1 had very different metal abundances: [Si/H] 
= -2.7 and -1.9. Toward Q1009+2956 Buries & Tytler 
(1998) found that two components also separated by a sim- 
ilar amount had abundances of [Si/H] = —2.7 and —2.4. 
Hence for these two QSOs the metal lines do not give a 
sufficient description of the H I and D I. Toward DLAs 
with log N H i> 20.3 cm" 2 , Wolfe & Prochaska (2000) find 
that different low ionization ions show similar velocity dis- 
tributions, which implies that both the ionizations and rel- 
ative abundances are similar, but does not constrain the 
absolute abundances. In those cases the differences were 
in different components, but for Q1243+3047 the veloc- 
ity offset would require structure within component 1, for 
which we have no other evidence. 

5.1.3. Velocity Agreement for D-2 

We suspect that D-2 is D for three reasons, but the case 
is weak, because D-2 is not well determined. First, the ab- 
sorption that we fit with D-2 is very narrow, narrower then 
typical H I lines. Second, the separation of D-2 from D-l, 
12 km s -1 , with a ler range of 8.5 - 14 km s _1 , matches the 
separation of the O I components, that is ~ 13.3 km s _1 . 
In the grid search we examined D separations from 4 - 
20 km s" 1 , and the D separation could have been differ- 
ent from that of the O I. However, this argument is weak, 
in part because the D-l and O-l are at slightly different 
velocities. Third, the fraction of the D in D-2 (0.12^1) 
matches the fraction of O in 0-2, 0.13, when we force 6(0- 
1) = 6(0-2). 

5.2. Line Width Agreement 

The widths of the D I, H I and O I lines indicate that 
the D-l absorption is D. If the D I, H I and O I arise 
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in the same homogeneous gas, and the velocity distri- 
bution in this gas is described by a single temperature, 
and a Gaussian distribution of turbulent velocities, then 
we can predict the 6-value of the D from the H and O 
b values. Following O'Meara et al. (2001), the intrin- 
sic b value of the lines, b int is given by b\ nt = b 2 emp + 
bturb> wnere the temperature term is b 2 emp = 2kT/m = 
166.41(km s _1 ) 2 (T/10 4 K)/mass (amu), where T is the 
gas temperature and bturb represents the bulk turbulent 
motions, and m is the mass in atomic mass units. The b 
values in this paper are all intrinsic values, since we have 
convolved the intrinsic line profiles by the instrumental 
broadening before we fit to the spectrum. 

The observed width of the D-l component agrees with 
the prediction from the widths of the 0-1 and H-l compo- 
nents, given in Table 3, and shown in Fig. 14. The straight 
line that fits b 2 against 1/m gives T = 1.1 ± 0.6 x 10 4 K 
and bturb = 5.8 ±0.6 km s -1 , both reasonable values, sim- 
ilar to those that we found in other D/H systems. This 
fit predicts 6(D-1)= 11.3 ± 1.8 km s _1 , 1.2a larger than 
the observed 6(D-1)= 9.2 ± 0.2 km s _1 . This agreement 
suggests that D-l is D and not H. However the evidence is 
weak because the 6(H-1) value and its error, from the co- 
variance matrix, are both poorly constrained, and because 
the velocity difference of the D-l and O-l components may 
mean that the model we used to predict 6(D-1) is too sim- 
plistic. 

When we accept that the D-l component is D, we can 
use D-l together with the H-l and O-l to improve our 
estimates of the gas properties. We find T = 0.55 ± 
0.04 x 10 4 K, which is low but reasonable for the large 
log Njji, and bturb = 6.3 ± 0.2 km s _1 . The data arc 
also fully consistent with this fit, which we show as the 
dashed line in Fig. 14. Compared to this line, the mea- 
surements give x 2 ~ 1-38 for one degree of freedom, where 
Prob(x? > 1.38) = 0.24. 

6. MEASUREMENT OF THE H COLUMN DENSITY 

In this section, we describe the measurement of the 
Hydrogen column density in the system that shows D, 
along with new methods that we developed to improve 
the accuracy and reliability of Nhi measurements. We first 
find the approximate log Nhi value (§6.1), then a detailed 
model that fits the spectrum (§6.2), and then we examine 
other models that give acceptable fits (§6.3). We find the 
log Nhi value and its error in §6.4, and we discuss these 
values in §6.5. 

Most information on the Nhi comes from the shape of 
the damped Lya line at 4285 A that we will call the DLA 
because it shows damping wings, although its log Nhi 
value is less than the usual definition for a DLA system. 
Compared to the methods we used to get log Ndi, for the 
log Nhi we use different ways to fit the continuum and 
Lya forest, to explore the parameter values and to assign 
an error to log Nhi- 

6.1. The Approximate log JVhi Value 

We can quickly establish that the DLA has log Nhi— 
19.7 cm -2 . If we fit the HIRES spectrum by a Lya line 
with a lower column density, there remains unexplained 
absorption which looks like damping wings. This is very 
clear when log Nhi< 19.5 cm -2 , and can still be seen 



at larger log Nhi- On the other hand, when we fit with 
log Nhi= 19.8 cm -2 wc absorb flux which is seen, which 
is un-physical. To allow this extra absorption, the QSO 
continuum, including any emission line flux, would have 
to bend upward, on either side of the damped line, near 
4288 A and 4293 A, which is best seen when we divide the 
spectrum by the line profile. We then know that 19.5 < 
log Nhi cm -2 < 19.8, where the larger limit is twice the 
smaller, a range that is 20 times the ler range on the 
log N D i (total). 

6.2. An Initial Model for the Spectrum Near the DLA 

We now describe how we made an initial model of the 
spectrum. In §6.3 we will use this model as a starting 
point when we explore alternative models. 

A model of the flux in the spectrum is determined by 
three main factors: the flux emitted by the QSO, the Lya 
forest and the DLA. We developed our own software to 
control these coupled factors. We can change the model 
manually, or automatically, optimizing over parameters to 
give models that have the smallest \ 2 difference from the 
spectrum. These tools are described in Appendices I and 
II. 

To construct our model, we first placed a DLA at 
z = 2.526 with log Nhi = 19.7 cm -2 . We then set a pre- 
liminary continuum level from 4200A to 4350A, such that 
the continuum plus DLA touched most of the peaks in 
the spectrum. This continuum was defined by 10 control 
points: six between 4320A and 4340A to define the top of 
the Lya emission line, two at A < 4250 A to define the 
continuum far from the emission line, and two between 
4260 A and 4325A nearest to the DLA. This continuum 
has enough freedom to take on a wide variety of shapes in 
the crucial region where the DLA is strongly absorbing. In 
addition, we add Lya forest absorbers where we see lines. 
We then made made small modifications to the contin- 
uum, the log Nhi of the DLA and the Lya forest, until we 
arrive at a model that appeared to accurately reproduce 
the spectrum. This model used log Nhi = 19.695 cm~ 2 . 
Finally, we used this model as the starting point for an au- 
tomatic optimization that returned log Nhi = 19.73 cm~ 2 . 
This model, shown in Figure 15, has a reasonable x , a 
smooth looking continuum, the Lya forest absorbers are 
not unusual, and there is little discernible structure in the 
residuals. 

The error on the log Nhi from the optimization covari- 
ance matrix, er(log Nhi) = 0.0002 cm -2 , is too small for 
at least two reasons: (1) the final covariance matrix has 
significant off-diagonal term involving the Nhi, and (2) 
the x 2 manifold may have multiple minima. We believe 
that errors derived from the covariance matrix will be too 
small for all similar models that contain many correlated 
parameters. We discuss this issue in §14.1. 

The initial model has two problems. First, we do not 
know the error on the log Nhi value, and second, we have 
no reason to believe that it is unique. 

6.3. Restarting the Optimization to Find the Nm Range 

We now try to determine the robustness and uniqueness 
of our fits by seeing if, for a variety of input values for the 
parameters, we return to our best value of log Nhi = 19.73 
cm -2 . A fit is robust if it can be found when wc start from 
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a wide variety of input parameters. A fit is unique if there 
are no other fits which have significantly lower % 2 . Al- 
though we can readily find local minima in the manifold 
of parameters, it is well known that there are no simple 
ways of showing that we have found the unique, or global 
minimum. On the contrary, we expect that there are other 
fits, which may have different log Nhi values, and similar 
or lower x 2 values. 

We used the optimizing code to searched for fits with low 
X 2 values. We performed thousands of optimizations, each 
beginning with different parameter values. We added ran- 
dom numbers, selected from normal distributions (NDs), 
to the parameter values of our initial model (§6.2). The 
ND for the continuum control points had a = 1 unit of 
flux. Since the continuum near the damped line is 12 
units of flux, these input continua usually have significant 
bumps and dips. The ND for the log Nhi, b, and z of each 
Lya forest line had a a equal to the prior measurement 
error. 

We allowed the optimizer to move the two continuum 
points between 4260 A and 4325 A (shown as squares in 
Figure 15) in both wavelength and flux. The optimizer 
also varied the flux, but not the wavelengths of all the 
other continuum points. Some examples of the starting 
continua are shown in Fig. 16. 

We found that the velocity structure of the H I has 
little effect on the log Nhi value. For all optimizations, 
we distributed the total log Nhi of the DLA between the 
two velocity components seen in the O I line, and we 
put 13% of the total log Nhi in the second component at 
v = 13.3 km s . This proportion is similar to that seen in 
O I and D I. However, we measured the same total log Nhi 
if all of the H I is in one component. Our fits to the DLA 
give log Nhi (total) that excludes the H I from components 
3, 4 and 5. We fit these other components separately with 
the results listed in Table 3. We did not fit separate com- 
ponents near 20 and 40 km s _1 because they are poorly 
constrained, none are needed near 40 km s -1 , and gas in 
this range should have log Nhi< 17 cm -2 . 

The total log Nhi was not allowed to vary in any of the 
optimizations. However, the input log Nhi was different 
in different optimization re-starts. 

We conducted thousands of re-starts for each chosen Nhi 
value. The \ 2 a t the starts were typically 10 5 , while at 
the end, depending on the Nhi value, 5700-6700, similar 
to the degrees of freedom. We found reasonable solutions 
with low x 2 values for a wide variety of Nhi values. In 
Figure 17 we show the x 2 values for these re-starts. 

Some of the final solutions found by the re-start pro- 
cedure had Lya forest absorbers that are not typical of 
those seen in the Lya forest of other QSOs (Pettini et al. 
1990; Hu et al. 1995; Kirkman & Tytler 1997). In par- 
ticular, many of the solutions with log Nhi < 19.68 cm~ 2 
had very wide Lya forest absorbers near to 4295A. We do 
not believe that these solutions indicate that there may 
be wide Lya forest absorption near to the LLS, but rather 
that the optimizer can, if it needs to, use Lya forest ab- 
sorption with large b values to compensate for low values 
of Nhi in the LLS. For this reason we have rejected, from 
the bottom panel of Figure 17, all fits that included Lya 
forest lines with b > 150 km s _1 . In Figure 18 we show 
the effect of one broad Lya forest line on a portion of the 



spectrum. 

6.4. The Best Estimate for log Nm and its Error 

We now discuss our best estimate for the log Nhi value 
and its error. 

We found that the re-starts with log Nhi= 19.73 cm~ 2 
consistently gave the lowest x 2 values. If our re-starts had 
fully explored the parameter space, we would have con- 
cluded log Nhi= 19.73 ± 0.005 cm~ 2 , where the error is 
the range of log Nhi values have fits with x 2 < Xmin + ^ 
the Sx 2 method that we used for the log Ndi in §4. For no 
other values of log Nhi, including 19.72 and 19.74 cm -2 
did we find fits with x 2 within 10 of the minimum, let 
alone 1. 

We do not use this method to estimate the log Nhi error 
because our optimization process is not efficiently finding 
fits with the lowest x 2 values. For each log Nhi value, 
Figure 17 shows a range of x 2 values of approximately 
200, and the lowest x 2 is often 10 lower than the second 
lowest. To define the error on log Nhi using Sx 2 = 1, we 
would require many fits with x 2 values within 1 of the min- 
imum, and we would also have to prove that we had ade- 
quately explored all relevant models and parameter ranges 
(§4.2.8). 

We determine the error on the log Nhi value by testing 
the hypothesis that there is at least one acceptable fit for 
a given log Nhi value. We accept an log Nhi value if the 
restarts include one or more models with an acceptable x 2 
value. We accept a x 2 if there is > 5% probability of a 
larger value when the hypothesis, that the data came from 
the model, is true. We fit 5048 pixels near Lya with 289 
parameters, leaving 4759 degrees of freedom and hence the 
maximum acceptable x 2 value is 4920.6. We increased this 
maximum acceptable x 2 by a factor of 1.23, to 6052.3, for 
reasons discussed in §6.5.4. 

We found acceptable fits for the re-starts for log Nhi = 
19.68 - 19.78 cm' 2 when we reject fits with Lya forest 
lines with b > 150 km s -1 , and 19.70 - 19.78 cm~ 2 when 
we reject b > 100 km s _1 . Therefore, we choose 19.68 - 
19.78 cm -2 as the range of acceptable log Nhi values. 

It is helpful to assign a probability distribution func- 
tion to summarize the range of likely log Nhi values. We 
made the following choices: a normal function, centered at 
19.73 cm~ 2 , with an 80% probability that the true log Nhi 
is in the range 19.68 - 19.78 cm~ 2 . We choose 80% be- 
cause we guess that there is of order a 20% chance that 
the true log Nhi is not in the range 19.68 - 19.78 cm~ 2 , 
and because we do not want the function to be too sharply 
centered on 19.73 cm~ 2 . 

Our best estimate for the column of the H I in compo- 
nents 1 and 2, that is associated with log Ndi in the same 
components, is then 

log N ffl (ioiaZ) = 19.73 ± 0.04 cm- 2 , (2) 

which is a la error of 9%, similar to the error on the 
log N D i. 

6.5. Discussion of Nhi Value and its Error 

We now discuss our log Nhi value and its error. We 
discussed a similar set of issues for the log Ndi value in 
§4.2. 
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6.5.1. The Coupling of the Continuum, Lya forest and 

the log TVhi 

The flux emitted by the QSO varies smoothly over scales 
of thousands of km s _1 . Its precise shape is unknown a 
priori and is hard to gage from spectra because of the 
absorption and emission lines. 

The Lya forest lines are narrow relative to the uncer- 
tainties in the shapes of the continuum and the DLA. We 
can identify and fit these lines individually, and we can 
distinguish their effects from the shape of the emitted flux 
and the DLA. It is unlikely that such lines would have 
the velocities and column densities which make absorp- 
tion that varies smoothly over hundreds or thousands of 
km s" 1 , however the optimizer will attempt to do this to 
accommodate errors in the log Nhi and the continuum. 

The damped Lya line has a very well determined shape 
that gives the log Nhi- The information on the Nhi comes 
from the whole profile of the damped Lya line. While 
the line is most conspicuous over the central 20 A, the 
damping wings extend much further, and continue to ab- 
sorb about 1% of the flux from approximately 4233 A to 
4340 A. This range covers approximately ±3700 km s _1 
in cither direction from the line center, or 3500 pixels in 
total, and extends over the peak of the QSO Lya emission 
line at 4331 A. 

In Fig. 19 we show the change in the flux corresponding 
to a change in log Nhi of 0.04 cm -2 , the la error that we 
quote. This is the derivative of the flux with respect to 
the column density, in different units. The change in the 
flux is largest ±5.4 A or 380 km s _1 from the line cen- 
ter, where the line absorbs 60% of the flux. The change 
in the flux decreases nearer to the line core because little 
flux remains there. The fractional change in the flux keeps 
increasing toward the line center, but we can not detect 
this. The information in the spectrum that gives log Nhi 
comes from the whole of the Lya line. The change in the 
flux is significant in the core where there is almost no flux, 
and also far from the core where we can integrate over a 
wide wavelength range. 

The error on the log Nhi value for the DLA is larger 
than would be thought from considering just the change 
in the spectrum produced by changing the DLA line alone, 
because the Lya forest and continuum can be adjusted to 
accommodate some of the effect of changing the log Nhi- 

Figure 19 shows the change in flux required of a fit to the 
spectrum when the log Nhi changes slightly. If we start 
with a model that gives a good fit to the spectrum and 
then increase log Nhi, we could maintain a good fit if we 
add bumps to the continuum approximately 5.4 A from 
the center of the DLA. For small change's in log Nhi the 
optimizer is able to avoid adding bumps, and instead it 
can find good fits using different smooth continua shapes, 
and different amounts of Lya forest absorption; less where 
the dips were expected and more on either side of them. 
The optimizer adjusts all parameters of the model so that 
the residuals appear uniformly distributed. This is one 
type of behavior that explains why the errors on log Nhi 
are greater than we expect by varying log Nhi alone. 

We now discuss in more detail the connection between 
the error on log Nhi and how we fit the Lya forest and the 
continuum. 



6.5.2. Fitting the Lya forest 

The models accommodate different log Nhi values by 
changing both the continuum and the Lya forest together. 
We must be careful to not allow the optimizer too much 
freedom to adjust the Lya forest parameters to simulate 
smoothly varying absorption, to accommodate an exces- 
sive range of log Nhi- For example, when we divide the 
spectrum by a DLA with log Nhi = 19.8 cm~ 2 we make 
a spike in the flux near 4293 A that is unlike any feature 
seen in QSO spectra. Hence, we should not allow the op- 
timizer to obtain acceptable fits for this log Nhi value by 
raising the continuum over a wide wavelength range and 
using smooth Lya forest absorption to remove this extra 
flux everywhere, except at the spike. This is reason that 
we require all of the Lya forest fe-values to be less than 
150 km s _1 . 

We may have over-estimated the error on the log Nhi if 
we have allowed the Lya forest lines to have more variety 
than is typical. The Lya forest can accommodate changes 
in the log Nhi in two ways: by making a few lines un- 
usually broad, and by using more lines than normal. We 
rejected the fits that used lines with b > 150 km s _1 , a 
value that is not well defined. If we instead reject fits with 
b < 100 km s _1 , we reduce the allowed range of log Nhi- 

Alternatively, we may have underestimated the error on 
log Nhi because we did not adequately explore the Lya 
forest absorption. We began all restarts using the param- 
eter values that deviated about the values of the initial 
fit. Perhaps the deviations that we gave to the Lya forest 
lines were too small to allow them to adequately fit the 
continum required for differing log Nhi- The deviations 
that we gave to the Lya forest parameters were uncorre- 
cted, so that increased absorption by one line would tend 
to cancel decreased absorption by its neighbors. We might 
have explored the Lya forest more thoroughly if we had 
used correlated deviations. 

We have not performed any quantitative checks of 
whether we have used an accurate Lya forest absorption 
model, but the excessive x 2 value of our model (§6.5.4) 
suggest that we may have systematically used too few com- 
ponents in our Lya forest model. We suspect that the 
effect this has Nhi is systematic, though we do not know 
in which direction. We also suspect the effect is small, in 
part because two of us developed independent models of 
the Lya forest absorption, both of which gave the same 
log Nhi = 19.73 cm~ 2 . It's hard to proceed quantitatively 
with this issue because it is extremely difficult to develop 
alternative models of the Lya forest absorption (number of 
lines, their approximate positions, etc) which arc equally 
compelling to the one we used. 

6.5.3. Continuum Shape 

Because the DLA absorbs over a very large wavelength 
range, we suspect that the largest class of systematic er- 
rors which might affect the Nhi are those where the flux 
emitted by the QSO can not be adequately fit by our con- 
tinuum model. This might occur, for example, if there are 
weak emission lines or high frequency errors in the relative 
flux calibration, in the vicinity of the DLA. 

The continuum shape can be effected by errors in the 
flux calibration. These errors are approximately 3% near 
the DLA, from comparison of HIRES spectra that we cal- 
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ibrated in different ways. We find the same log Nhi= 
19.73 cm -2 gives the lowest \ 2 values when we fit the DLA 
in HIRES spectra that we calibrated in different ways. 

The models described above all used only two contin- 
uum control points that were free to move in both flux 
and wavelength near the DLA. We also performed re-starts 
with three such continuum control points in the region be- 
tween 4260 A and 4325 A. The results of these re-starts 
are shown in Figure 20. The top panel contains a number 
of fits that are unreasonable, having either very unusual 
Lya forest absorption (like some fits in the top panel of 
Fig. 17), or very bumpy continua that we could not make 
with only two continuum points. 

We filter the results using 3 continuum control points 
to exclude the unreasonable fits, leaving the results shown 
in the bottom panel of Figure 20. As with the two con- 
tinuum point re-starts, we first removed all fits containing 
Lya forest absorbers with b > 150 km s _1 . In addition, 
we removed models having continuum bumps larger than 
we expect to be present due to cither weak emission lines, 
broad shallow absorption or errors from the relative flux 
calibration. 

To help us decide what size of bumps and dips were 
reasonable in QSO continua, we fitted the Lya forest of 
three other QSOs in a similar manner to Q1243+3047. 
Our HIRES spectra of these QSOs had similar, or in two 
cases a bit lower S/N. Near the position of the DLA in 
Q1243+3047, approximately 3000 km s _1 from the peak 
of the Lya emission line, we did not require any bumps or 
dips in the continuum that exceeded approximately 2% in 
amplitude, similar to the error in the relative flux calibra- 
tion for Q1243+3047. 

The filter that we used to reject bumpy continua was 
guided by this result, but it was imperfect in its treatment 
of the Lya emission line. Examples of continua that it 
ejected are shown in Figure 21, and those that it accepted 
in Figure 22. 

We find that the results from the models with two- and 
three-control point in the continuum are similar. This 
suggests that we have explored an adequate range of con- 
tinuum shapes to accommodate the emitted flux, weak 
emission lines, broad shallow absorption and errors in the 
relative flux calibration. 

The first thing that we note is that the best model with 
three-continuum points has log Nhi = 19.73 cm -2 , the 
same value as with two control points. Second, the range of 
solutions is also approximately the same; log Nhi = 19.71 
to just above log Nhi = 19.78, which was unfortunately 
the highest log Nhi we explored with three points. 

The third finding is that most of the filtered 3-point con- 
tinua with acceptable x 2 had shapes similar to those from 
the 2-point fits. In Figure 23 we compare the continua 
from the two and three-point continua, for the models with 
that gave the lowest x 2 values for log Nhi= 19.69 cm~ 2 , 
the value ler below our estimate for log Nhi- They dif- 
fer by at most < 1% in flux, which is less than the 
error from the flux calibration. In Figure 24 we show 
the same for log Nhi= 19.783 cm~ 2 , our best estimate 
for log Nhi- The differences are at most ~ 0.3%. For 
log Nhi= 19.783 cm -2 , shown in Figure 25, the differences 
are < 3%. 

We use the 2-point rather than the 3-point continua to 



define the log Nhi because we have explored the former in 
more detail, and the latter give a large number of fits with 
unacceptably bumpy continua, that we can not readily fil- 
ter. 

Although the 2 and 3-point continua give very similar 
results, we are not certain that we have explored an appro- 
priate range of continua. We might have allowed either too 
much or too little freedom to represent the continuum, and 
we do not know how such errors might effect the log Nhi 
value and its error. 

6.5.4. The Largest Acceptable x 2 Value 

As we stated in §6.4, we increased the maximum accept- 
able x 2 by a factor of 1.23. We did this scaling because we 
found that when we fit a similar region of the Lya forest 
without unusually strong lines, we obtained a \ 2 value of 
6111 over 5290 pixels using 336 parameters. This is x 2 per 
degree of freedom is much larger than we would expect. 

One possible explanation is that we are under-fitting the 
Lya forest. This is reasonable because we fit lines only 
where they are clearly required. We did not fit lines that 
changed the residuals by less that approximately ler, and 
hence we have underfit the Lya forest because the need 
for such lines would be apparent where the S/N higher. 
We did not fit such lines because the spectrum gives little 
guidance in how to obtain a lower x 2 ■ The residuals, like 
those shown in Figure 15, do not show the need for distinct 
lines. 

A less likely explanation is that the error array associ- 
ated with the spectrum is systematically too small. When 
we fit 5882 pixels to the red of the Lya emission lines, 
using just 32 parameters, we found a x 2 value of 5264, or 
0.90 per degree of freedom, which suggests that the errors 
are on average too large. 

Several factors will determine the \ 2 value we obtain 
when we fit the Lya forest. The x 2 w iU decrease when 
we add lines. As the S/N increases, we need more lines to 
reach a given x 2 per pixel, and the higher S/N helps show 
that these lines are needed. The S/N decreases system- 
atically with decreasing wavelength in our HIRES spec- 
trum, and hence we do not know if the 1.23 factor applies 
at all wavelengths. In addition, the Lya forest evolves, 
with more absorption at higher rcdshifts, until we near 
the QSO, where the proximity effect reverses the trend, 
further complicating the issue. 

Although we can not readily determine the accuracy 
of the correction factor of 1.23, Fig. 17 shows that the 
log Nhi range is rather insensitive to the precise maximum 
acceptable x 2 value. 

6.5.5. Speculation on the Errors Associated with the x 2 

Range 

The portion of the error on the log Nhi that comes 
from only the quality of the fit and the S/N of the spec- 
trum, might be up to an order of magnitude smaller than 
0.04 cm -2 . We speculate this because the 0.04 cm -2 error 
corresponds to a range of x 2 of 250, while we could have 
used a range of x 2 °f 1 had we adequately explored both 
the parameter and model space. Because we have not done 
the required exploration, we do not know the shape of the 
Xmin as a function of log Nhi to the required accuracy. 
If we make the guess that the shape remained similar to 
that suggested in Fig. 17, then the range of log Nhi corre- 
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sponding to Sx 2 = 1 would be much less than 0.04 cm" 2 . 
Of course, the central value could end up anywhere in the 
range indicated by Equation 2. 

7. ABUNDANCES AND IONIZATION OF THE HEAVY 
ELEMENTS IN THE ABSORBING GAS 

In this section we discuss the metal abundance, ioniza- 
tion, and physical conditions in the components (1 and 2) 
where we measure D/H. In §7.3 we discuss the components 
at -40 and 100 km s" 1 . 

We modelled the level of photoionization using 
CLOUDY version 94.00, developed by G. Ferland (Fer- 
land, 1991). We used the solar Oxygen abundance log 
O/H = —3.31 from Allendc, Prieto, Lambert & Asplund 
(2001). For other elements we used the CLOUDY de- 
fined solar abundance ratios. We assumed a plane-parallel 
geometry and we approximated an isotropic background 
by placing a point source at a very large distance. We 
used the Haardt-Madau (Haardt & Madau, 1996) ioniz- 
ing spectrum at z = 2.526. The ionization is given by 
the parameter logt/, where U = <p/cnjj = J912/ ^hcriH, 
(j> (cm -2 s" 1 ) is the surface flux of ionizing photons, and 
J912 is the intensity of the incident radiation at 1 Ryd- 
berg. We set J912 = 10 -21 ergs cm -2 s" 1 Hz -1 sr" 1 , a 
typical, if slightly high value for the intensity of the inter- 
galactic radiation field (Scott et al. 2000; Hui et al. 2002; 
but see Prochaska 1999). The gas density, n# cm -3 , is 
then related to the ionization by log U = —4.34 — log nn 
(cm -3 ). 

7.1. Metal Abundance in Components 1 & 2 

Photoionization models show that the gas in compo- 
nents 1 & 2 has relatively low ionization and hence we 
used O I/H I to find [O/H] ~ -2.79. 

We ran a set of CLOUDY models for log N H i= 
19.73 cm -2 and metal abundance [X/H] = —2.77, similar 
to the value we expected to find based upon the observed 

I/H I ratio. We considered all the gas in components 

1 and 2 together because there are large errors on the 
fraction of each ion in component 2. Figure 26 shows the 
predicted column densities for ions of interest as a function 
of gas density. 

The absence of strong C IV and Si IV shows that the 
density is high enough that the metal abundance can be 
taken from the O I/H I ratio. The components that we 
fit to C IV and Si IV nearest to v = km s" 1 have 
b ~ 38 km s _1 , very different from the b ~ 7 km s" 1 
that we expect for components 1 & 2. In Table 3 we list 
upper limits on the column densities that could be asso- 
ciated with components 1 & 2, obtained when we fit lines 
with fixed 6=7 km s" 1 and v = km s" 1 . The spec- 
trum allows these column densities, but not much larger 
values. We also list upper limits that we found for Si III 
and Al II. These limits all require relatively high densities: 
\ogn H (cm" 3 ) > -2.3 (Al II), -1.7 (Si IV), -1.6 (Si III) 
and -1.5 (C IV). 

We obtain the [O/H] from the observed O I/H I, making 
a slight adjustment for ionization. The O I column density 
in components 1 & 2, 13.63 ±0.02 cm -2 , equals the model 
prediction for [O/H] = —2.77 and log Uh = —2.24 cm -3 . 
Since the density is relatively high, the [O /H] is insensitive 
to the precise density. The other ions prefer a higher den- 



sity of lognH — — 1.5 cm" 3 , where the predicted O I col- 
umn density is 0.02 cm" 2 higher than we observe. Hence 
we conclude [O/H] = —2.79 ±0.05, which includes the er- 
ror on log Nhi- Had we chosen a much higher density of 
\ognn = +1 cm" 3 , the [O/H] would decrease by only 0.01. 
We checked that that the [O /H] value is insensitive to the 
shape of the ionizing spectrum. 

7.2. Ionization, Size and Mass of Components 1 & 2 

Overall, the column densities of the ions indicate 
log nn — —1.5 cm -3 , or \ogU ~ —2.84. Although a lower 
density of ~ —2.1 cm" 3 is preferred by Si II and is allowed 
by Al II, the C IV and Si IV rule this out, and the C II 
prefers a higher density, log nn ~ —0.44 cm" 3 . We used 
the column densities in the components at the velocities 
that we footnote in Table 3. The fit is improved if we 
add C II at v = —5.6 km s _1 , or if C is slightly under- 
abundant and Si is over- abundant, but none of these dif- 
ferences are significant because the column densities are 
not well known. 

The level of ionization indicated by the ions implies that 
the gas is mostly ionized. Figure 27 shows the neutral frac- 
tion is of order 21% for logU = -2.84. 

The log hh value implies an absorbing region of typical 
size and mass. The bottom half of the Fig. 27 shows the 
size for a given ionization, for either constant nn or con- 
stant Jgi2- For a given log U, nu is proportional to J912 by 
definition, the size is proportional to njj 1 , and the mass to 
nj{ 2 . For an ionization of logt/ = —2.84 and J912 = 10" 21 
ergs cm" 2 s" 1 Hz" 1 sr" 1 , the size is 2700 pc along the line 
of sight, and the mass of H in a sphere with this diameter, 
and the density logra# = —1.5 cm" 3 is 8.0 x 10 6 M Q . 

Both the size and mass of the absorbing region are un- 
certain by at least an order of magnitude because J912 
and especially the gas density are uncertain. Consider 
three types of change we can make to lower the sizes and 
masses: we can increase the density, we can decrease the 
ionization, or we can increase the density and the J912 in 
proportion to maintain the ionization, which we prefer. 
For example, if we increase both the J912 and the nu by 
10 times, the ionization is unchanged, the size drops by 10 
times, and mass by 100 times. 

If we consider low ionization solutions, we can have high 
gas densities that give small sizes and masses. We view 
this alternative as unlikely, because at lower ionization we 
expect less Si II than we see, and less Al II than allowed by 
the upper limit. This is possible if most of the absorption 
from these ions is from different gas. However, this would 
also require C II to come from different gas, yet Fig. 4 
suggests that much of the O I, C II and Si II do some from 
the same gas, and the ionization model is consistent with 
this. But if we do consider lower ionization solutions and 
we keep the density constant at log?i# = —1.5 cm" 3 , the 
size comes down by at most a factor of 5, and the mass by 
a factor of 125. This is because the gas becomes neutral, 
at which point the size depends only on the density. 

The properties of the gas that shows D towards 
Q1243+3047 are typical. Table 4 summarizes the inferred 
total gas column densities, sizes and gas masses of six D/H 
absorbers. For consistency we scaled the gas densities, and 
hence sizes and masses so that the observed ionization was 
obtained for J912 = 10" 21 ergs cm -2 s _1 Hz" 1 sr" 1 . We 
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see that the total gas column densities, log H, may all be 
very similar. The sizes are of order 1 kpc. but with a huge 
range covering over a factor of 1000, while the masses cover 
an even larger range. 

Finally, CLOUDY provides estimates for the gas tem- 
perature that we show in Fig. 28. For \ogU ~ —2.84 we 
expect an equilibrium gas temperature near 10,000 K. At 
this temperature we expect b(D-l) > 9.1 km s , similar 
to the observed value. We saw in §5.2 that the widths 
of the O-l, D-l, and H-l components together indicate 
T = 0.55±0.04 x 10 4 K, which is lower, perhaps because 
the gas is, after all, neutral or the ionizing spectrum is 
softer than the one we used. 

7.3. Abundances in other components 

Photoionization calculations suggest, with large uncer- 
tainty, that components 3 near v ~ —40, and 4 near 
v ~ +100 km s _1 , have different metal abundances from 
each other and from components 1 & 2. 

We performed additional photoionization calculations 
for a metal abundance [X/H] = —1.5 and njj = 0.001 
cm" 3 and log Nhi = 16.00 cm~ 2 (component 3) and 
log Nhi= 16.50 cm~ 2 (component 4). The nn was de- 
termined using the relation between column density and 
physical density as given by Hellsten et al. (1998). Figure 
29 illustrates the results. 

For component 3 we estimate [C/H] = —0.5 to —1.5, 
from Si II/C II and the lack of Si IV and C IV. The ion- 
ization is in the range — 5 < \oglf < —2.75, and the Si 
is not enhanced. The equilibrium temperature requires 
13.36 <b H < 17.50 km s"\ when we assume a turbulent 
velocity width of 3 km s -1 , as implied by the widths of the 
Si II and C II absorption. This range is consistent with 
the grid search value b(H — 3) = 16.3 ± 0.7 km s _1 (§4.1). 

For component 4 we estimate [Si/H] ~ [C/H] ~ —1.9 
with a large uncertainty, because different ions suggest dif- 
ferent ionization. If Si is enhanced by 0.3 dex, C Il/Si II 
suggests \ogU ~ -1.8 and [Si/H] ~ [C/H] ~ -1.7. How- 
ever, C IV/C II and C IV/Si IV indicate higher ionization, 
\ogU ~ -2.6, and [Si/H] ~ -2.1. 

Three of the components of the absorption system at 
z ~ 2.526 appear to have very different abundances: —2.79 
for the gas in components 1 & 2, —2 for component 4 and 
— 1 for component 3. We note that the component with 
the most H I and the lowest ionization also has the lowest 
metal abundance. 

8. BEST FIT VALUES FOR Q1243+3047 

We have made the following estimates for the parame- 
ters that describe the gas in which we measure D/H: 

• The velocity structure of the gas has two main com- 
ponents that we see in the asymmetry of the O I 
line. 

• Our grid search of models that fit the D lines showed 
that the components are separated by 12 km s _1 , 
with a ler range of 8.5 - 14 km s _1 , consistent with 
the asymmetry of the O I line. 

• The main D component, D-l, is at z = 2.525659 and 
has b = 9.2 ±0.2 km s _1 . 



• The wavelength, 6-value and column density show 
that D-l is D absorption. We believe that the D-2 
component is D because it is very narrow, the sep- 
aration of the D-l and D-2 components is the same 
as that of the O-l and 0-2 components, and the 
fraction of the D in D-l is similar to the fraction of 
the O in O-l (§5). However, these arguments are 
weak because D-2 is not well defined. 

• The grid search also indicated that component D-2, 
at the larger redshift, contains 12% of the total D 
column density, with a range of 7% to 28% depend- 
ing on how much of the absorption near component 
2 is H in another component, H-3. 

• Total logN ra = 15.1131^26 cm" 2 , (i.e. +10%, 
—6%) where the error range includes fits with x 2 < 
Xmin + 1- The errors are larger for higher values of 
log Ndi because of the second component can make 
a significant contribution. 

• Total log N ffl = 19.73 ± 0.04 cm" 2 (9%), where the 
error is a Gaussian distribution centered on the 
log Nhi value that consistently gave the best fits. 
We set the width of the Gaussian to give an 80% 
chance that the true log Nhi lies in the range 19.68 
- 19.78 cm" 2 where we found acceptable fits to the 
Lya line, the continuum and Lya forest. 

. log (D/H) = -4.6171S (+14%, -10%). 

• D/H = 2.42±°; 3 JxlO~ 5 cm" 2 . 

• [O/H] = -2.79 + 0.05. 

• The absorber is probably mostly ionized, with of 
order 21% of H atoms neutral (§7). 

9. THE PRIMORDIAL d/h RATIO 

In this section we compare D/H measurements from dif- 
ferent QSOs, we discuss the dispersion in these value, and 
whether D/H might correlate with metal abundance or 
Nhi, and we give our estimate for the primordial D/H 
value. 

9.1. The Weighted Mean D/H from Five QSOs and the 
Dispersion of the Values 

We find the weighted mean of the D /H values from five 
QSOs and we show that the individual values show more 
dispersion than we expect. In Tables 5 & 6 we list all re- 
ported D/H measurements that remain viable. We have 
previously measured D/H in three QSOs (Tytler, Fan & 
Buries 1996; Tytler & Buries 1997; Buries & Tytler 1998a; 
Buries & Tytler 1998b; O'Mcara et al. 2001), and placed 
a strong upper limit on D/H in a fourth (Kirkman et al. 
2000). We will also use the D/H measurement by Pettini & 
Bowen (2001) towards Q2206-199, although this measure- 
ment is less secure because the HST spectra are of much 
lower S/N and resolution than those from the ground. We 
discuss Q2206-199 and QSOs that we will not use in Ap- 
pendix III (§15). 

The weighted mean of the first five log D /H values from 
Table 6 is, 

log D/H = -4.556, (3) 



15 



where the weights we use are the la errors on the quantity 
Yi = log{D / H)i. We use log values because they were used 
to find all but one of the individual D/H values and errors. 
We obtain a slightly smaller mean D /H if we instead work 
with the linear D/H values. 

The D/H measurements towards the five QSOs are more 
dispersed than we expect. In O'Meara et al. (2001) we 
noted that the dispersion in the first three measurements 
was larger than expected, with a 3% chance of a larger x 2 
value. We interpreted this to mean that we had underes- 
timated one or more of those errors. With the addition 
Q 1243+3047 the dispersion of log D/H remains approx- 
imately 0.10, but adding the low D/H from Q2206-199 
increased the dispersion to 0.14, the x 2 value for all five 
measurements increases to 12.35 for 4 degrees of freedom, 
and the probability that we would have obtained a larger 
X 2 value by chance drops to 1.5%. 

9.2. Factors that Determine Measurement Errors on 
D/H Values 

The accuracy and reliability of a D /H value will depend 
on many factors, that fall under three main headings: 

• The details of the absorption system: the Nhi value, 
the number of velocity components, their velocities, 
column densities and b- values, the chance placement 
of Lya forest lines and other contamination, the red- 
shift and the flatness of the continuum near the key 
lines. 

• The quality of the spectrum: the spectral resolu- 
tion, the S/N, the ions observed, and the accuracy 
of the wavelengths and relative flux calibration. 

• The adequacy of the model: Have all possible veloc- 
ity structures been explored? Have alternative line 
identifications been considered? Where contaminat- 
ing and blended lines fit? Were the hidden compo- 
nents explored? Is the continuum over- or under-fit? 
Have the continuum, Lya forest and the lines that 
give D/H been fit simultaneously? 

The dominant error with any particular D/H value 
might depend on any of these factors, or some combina- 
tion of them. Many of these factors are different between 
the D/H values, and hence each value can have a differ- 
ent dominant error. We do not expect the measurement 
errors, or any uncorrected systematic errors to correlate 
with a single parameter, such as the z a b s , Nhi or metal 
abundance. Although unlikely, we do not know enough to 
rule this out. 

The errors that are associated with any D/H value could 
be too small for any of the factors in the list above, such 
as: poor flux calibration, especially near the Lyman lim- 
its, inadequate fitting of the Lya forest, and inadequate 
modelling of the continuum and velocity structure, and 
especially, inadequate exploration of all possible models. 
The error on a D/H value could be larger than that de- 
duced from the % 2 values for any of these reasons. 

9.3. The Dispersion in the D/H Values May Come from 
Measurement Errors 



We suspect that the dispersion in the D /H values arises 
from measurement errors and is not real. If the mea- 
surement errors have been underestimated for at least one 
QSO, then we can explain the excess x 2 value. 

The dispersion of the D /H values is not much larger than 
we expect. In Table 6 we list Xi — (Yi — mean)/cr(Yj), the 
deviation of each measurement from the weighted mean, in 
units of the individual measurement errors. The D/H from 
Q1009+2956 is 1.95cr above the weighted mean, while the 
D/H from Q2206-199 is 2.17cr below. A X 2 < 9.5 would 
have been expected if all five log D/H values were consis- 
tent with the weighted mean, since Prob(%4 > 9.5) = 0.05. 
We would obtain X 2 < 9.5 if either Q1009+2956 or Q2206- 
199 were within la of the mean, or if the measurement 
errors on all five QSOs were increased by 1.14, both small 
changes. However, to obtain a typical x 2 = 3.36, where 
Prob(x4 > 3.36) = 0.5, 3 or more D/H values, or their 
errors, would need to change, or the measurement errors 
on all five QSOs would need to increase by a factor of 1.92. 
These arc a large, but still credible changes. 

The methods that we introduced in this paper have al- 
lowed us to explore some of the issues that effect a D/H 
value more thoroughly that in past work. We sec many 
ways in which the errors might have been underestimated, 
both here, and in past work. In §4.2 and §6.5 we discuss 
factors that might change the la error that we give for 
the D/H value from Q1243+3047 by a factor of two. We 
further believe that the D/H value for Q1243+3047 has 
benefited from the most thorough exploration of these er- 
ror related issues. We would not be surprised if the errors 
on some of the other D /H values were too small by a factor 
of two. 

We are now more confident that measurement errors 
explain the dispersion in D/H values than we were in 
O'Meara et al. Our work on Q 1243+3047 better illus- 
trates the many ways in which errors can be underesti- 
mated. We see that the three D/H values that we sus- 
pect are the least reliable are also the values farthest 
from the mean. For PKS 1937-1009 and Q1009+2956 
some of the data reduction and analysis methods were less 
reliable than those that we used for HS0105+1619 and 
Q1243+3047, and for Q2206-199 there is much less data 
available. 

9.4. How D/H Depend on Metal Abundance: D/H 
Chemical Evolution 

The mean D/H from QSOs is similar to the primordial 
D/H value that we would predict from the D/H and metal 
abundance in the local interstellar medium (LISM) using 
standard Galactic chemical evolution. However, chemical 
evolution can not account for the dispersion in the D/H 
values from QSOs. 

With the two latest D/H measurements, there is no 
longer a hint of a correlation between D/H and metal 
abundance in the QSO absorbers that we noted in 
O'Meara et al. (2001). Prantzos & Ishimaru (2001) showed 
that standard chemical evolution could not reproduce the 
correlation, while Fields et al. (2001) discussed an uncon- 
ventional scenario that could. In Fig. 30 we show D/H 
values against metal abundance, mostly measured with Si 
when the ionization is high, and O otherwise (Table 6). 

With the elimination of the correlation, Galactic chemi- 
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cal evolution more clearly supports the idea that our D /H 
measurement towards Q1243+3047, and the mean D/H 
from five QSOs are very close to the primordial D/H value. 
We include two curves in Fig. 30 that show the expected 
decrease of D/H against metal abundance in a simple 
closed box model, without in-fall, and with the instan- 
taneous recycling approximation (Tinsley 1974; Tinsley 
1980; Ostriker & Tinsley 1975). We normalized the two 
curves in different ways. The solid curve uses the mean 
D/H from Eqn. 3 as the primordial D/H and predicts 
D/H as a function of metal abundances, including that in 
the LISM. The dotted curve is normalized to give the D/H 
abundance in the LISM (Moos et al. 2002; Oliveira et al. 
2002) and it then predicts the primordial D/H abundance, 
(D/H) p . To draw these curves, we use the equation 

D/H = (D/H) p x cxp (~-JLSj , (4) 

and the values for two parameters from Prantzos & Ishi- 
maru(2001); the returned mass fraction R = 0.31, and 
the yield y = 0.6, where Z is the metallicity, or metal 
abundance in solar units, on a linear scale. These pa- 
rameters were derived from a Salpeter IMF and the yields 
of Woosley & Weaver (1995). The choice of parameters 
and D/H chemical evolution are discussed in Steigman 
& Tosi(1992), Vangioni-Flam, Olive & Prantzos (1994), 
Galli & Palla (1995), Prantzos (1996) and Prantzos & Silk 
(1998). 

Both curves connect the D/H values towards 
Q1243+3047 and HS 0105+1619 to the D/H in the LISM, 
over three orders of magnitudes in metallicity. The curves 
show that the decline in D/H, as stars eject gas that lacks 
D, is insignificant when metal abundances are low; e.g. 
[O/H] < —1. The simple closed box model (Eqn. 4) 
predicts that by the metal abundance of Q1243+3047, 
(D/H) = 0.9987(D/H) p , and for HS 0105+1619, 0.986. 
We have not applied these corrections because they are 
much smaller than our measurement errors. 

We predict the primordial D/H when we normalize the 
simple model to give the D/H in the LISM. Ostriker & 
Tinsley (1975) predicted the primordial D/H would be 1.5 
- 2 times that in the LISM. With the modern param- 
eters given above, the same model predicts 1.62 times, 
again without using any D/H measurements. Using the 
D/H and Z measurements in the LISM from FUSE and 
HST spectra, D/H LISM = 1.52 ± 0.07x 10~ 5 , and [O/H] 
= -0.189 or Z = 0.647 (Oliveira et al. 2002, and using 
the solar log O/H = —3.310 from Allcndc Pricto, Lambert 
& Asplund 2001), the predicted primordial abundance is 
D/H = 2.47+ 0.13xl0 -5 , very similar to the value from 
Q1243+3047 and from Eqn. 3: 2.78xl0~ 5 . The errors 
on the prediction are substantial. When we allow the 
two parameters to simultaneously take values that max- 
imize the change in D/H, from 0.26 < R < 0.36 and 
0.5Zq < yield < 0.9Zq (Prantzos, private communica- 
tion), we find D/H = 1.94 - 3.16xl0 -5 . The range is 
further increased when we consider in-fall of gas to the 
Galactic disk (Lubowich et al. 2000) and dispersion of 
D/H in the LISM (Vidal-Madjar et al. 1998; Sonneborn 
et al. 2000). Hence, although the D/H values from both 
Q1243+3047 and HS 0105+1619 arc closest to the pre- 
dicted D/H, we can not use the chemical evolution model 
to rule out primordial D/H values that are suggested by 



the measurement to the other QSOs. 

9.5. Does D/H Depend on N m ? 

The D/H values appear to decline with increasing 
log Nhi- We do not believe that this correlation is real, 
because we suspect that the dispersion in the D/H val- 
ues is not real (§9.3). We will discuss potential systematic 
errors, and then potential cosmological and astrophysical 
origins for a correlation. 

In Figure 31 we see a clear trend of declining D/H with 
Nhi- This trend was apparent with just the three D/H val- 
ues discussed by O'Meara et al. (2001) and it was accen- 
tuated by Q2206-199 from Pettini & Bowen (2001). The 
trend rests upon the relatively high D/H for the two LLS 
(PKS 1937-1009 and Q1009+2956), and the relatively low 
D/H for Q2206-199 and Q1243+3047 that we will call 
DLAs. 

9.5.1. Systematic Measurement Errors that Depend on 

log N m 

We are not aware of any systematic error in the measure- 
ments of PKS 1937-1009, Q1009+2956, HS 0105+1619 
and Q1243+3047 that could readily account all of the 
trend with Nhi- Rather, we expect that potential errors 
are complex, and specific to each spectrum (§9.2). How- 
ever, we can think of three types of systematic error that 
might depend on Nhi- 

One possibility is that we have included absorption by 
H in the D measurements, making D/H too large. Hydro- 
gen absorption is more likely to appear like D in the LLS 
compared to the DLAs for three reasons. First, the LLS 
have weaker D lines, and hence small amounts of H can 
have a larger effect than with the DLAs. Second, in the 
LLS we see only the D Lya line, and hence we have less 
information on the D velocity structure than we have for 
the DLAs. Third, the b values for the H and D lines in the 
LLS are larger than those in the systems with higher Nhi 
values, which makes the D lines less different from com- 
mon H lines. A strong argument against this possibility is 
that the D lines in both PKS 1937-1009 and Q1009+2956 
had 6-values that were smaller than Lya forest lines, and 
that agreed with predictions from their H and metal lines 
that we collect in Table 5. 

A second possibility is that we have underestimated the 
Nhi for the two LLS. The method of measuring Nhi is dif- 
ferent for the LLS than for the systems with higher Nhi, 
hence this could be a systematic effect. 

A third possibility is that we have missed some of the 
components of D, making D/H too low in some QSOs 
(§4.2.6). D components at positive velocities can hide in 
the parts of the spectrum where the main H absorbs. This 
can happen in any system where the v distribution of the 
H I is poorly known, and the H absorption extends over 
nearly 82 km s _1 , enough that the H with large positive 
velocities has D that lies among the absorption of the H 
with negative velocities. 

We do not know whether missing components are more 
likely at low or high log Nhi- It would be header for 
this effect to go undetected in systems with the low- 
est column densities and unsaturated Ly series lines, like 
Q1009+2956. We might expect this effect to be worst in 
systems with the highest log Nhi because the H lines are 
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then widest. However two other factors make it easier to 
find components of the systems with higher log Nhi val- 
ues: the components then have enough gas to show metal 
lines, and the b values for the H I are smaller than for the 
LLS. 

9.5.2. Both LLS and DLAs can give Reliable D/H 

In general, both DLAs and LLS can give accurate D/H 
values, and both can have significant advantages. On the 
other hand, both LLS and DLAs can, in specific cases, give 
unreliable D/H values. If we had many more D/H mea- 
surements, we could determine observationally the fraction 
of DLAs and LLS that had given errors, but we can not 
estimate this today. 

The DLAs have the advantages that we see several D 
lines, the b values can be much lower than in the LLS or 
the Lya forest, the D lines are strong enough that they 
are less effected by the Lya forest alone, and we can see 
many metal lines that can help us understand the velocity 
structure. 

The LLS with the lowest log Nhi also have major ad- 
vantages. They have unsaturated H lines that give the 
velocity distribution of the H. We can also measure the 
log Nhi independently the drop in flux across the Lyman 
break, as we did for PKS 1937-1009 (Buries & Tytler 1997) 
and Q1009+2956 (Buries & Tytler 1998b). These two ad- 
vantages make the log Nhi much easier to measure than in 
other systems. We can also directly compare the velocity 
distributions of the H and D, proving that the absorption 
is D, and guaranteeing that we are measuring the D and 
H in the same gas, and not missing components of either 
the D or the H. These advantages are offset by the greater 
similarity between the D and H lines. 

9.5.3. Cosmological and Astrophysical Reasons why D/H 
Might Depend on log Nm 

If D/H were correlated with log Nhi, the explanation 
might involve inhomogeneous BBN, or the creation, re- 
moval or destruction of D. No such plausible mechanisms 
are known. 

If we ignore measurement errors, the correlation in Fig. 
31 suggests that the D/H ratio has a range of approxi- 
mately a factor of 2.4. 

Although the correlation is seen against log Nhi, this 
might not be the most readily interpreted parameter. In 
Table 4 we list approximate, order of magnitude estimates 
for the physical conditions where we measure D/H. The 
log Nhi values correlate with the ionization of the gas, 
and with the approximate density and size and mass of 
the absorbing region. The highest log Nhi values (low 
D/H values) tend to correspond to low ionization, higher 
gas density and smaller, less massive absorbing regions. 
Q 1243+3047 stands apart from the trend because it has 
both a high log Nhi and it is ionized, giving a large size 
and mass, perhaps the largest of any of the D/H regions. 
Most of the masses are large, 10 5 - 10 9 M Q , which makes 
it hard to make or destroy the D in local events, unless 
they effect a large fraction of the higher density portions 
of the universe. 

In SBBN the range in D /H values would correspond to a 
range of approximately 2 in the rj value. This range might 
arrive is the universe is has inhomogeneous r\ on large 



scales at the time of BBN (Kainulaincn, Kurki-Suonio, & 
Sihvola 1999). The scale of the inhomogeneity would need 
to be > 1 kpc, and perhaps, for the log Nhi values of the 
D/H absorbers, > 100 kpc (Mike Norman, private commu- 
nication), to avoid mixing before the time of observation. 
The scale is also but is limited to < 1 Mpc by the near 
isotropy of the CMB (Jedamzik & Fuller 1995; Copi et al. 
1998; Jedamzik 2002). The 1 Mpc scale corresponds to a 
baryon mass of 5.9 x 10 9 M Q that is larger than the typical 
masses of baryon in the gas showing D/H. Inhomogeneity 
can produce varying D/H, but the power spectrum of the 
fluctuations should be cut off at both small and large scales 
to prevent the overproduction of 4 He and 7 Li, and to avoid 
CMB constraints. The lack of variation of the 7 Li/ H in 
stars in the halo of our galaxy is a weak argument against 
inhomogeneities on these scales, since we expect that the 
gas that we see in QSO absorption systems is the type of 
gas that makes the halo stars. The argument is weak be- 
cause the 7 Li in the halo stars may not be the primordial 
value (§10.1). 

Dolgov & Pagcl (1999) present a different scheme (dis- 
cussed by Kurki-Suonio 2000; Whitmirc & Schcrrer 2000 
and Dolgov 2002) that is inhomogeneous on scales of 100 - 
1000 Mpc, small enough to give variation in D/H to QSOs, 
but large enough to avoid large variations in existing 4 He 
and 7 Li measurements that are all in local objects. They 
employ inhomogeneities of the different neutrino flavors 
that add up to give a constant total energy density, thus 
avoiding the CMB constraint. 

We can speculate that significant D has been destroyed 
in those parts of the universe that, by the time of ob- 
servation, had the largest H I column densities. These are 
regions with larger over-densities relative to the mean den- 
sity. Fields et al. (2001) discuss a highly unconventional 
chemical evolution models and find three constraints on 
the conditions required to destroy significant D while keep- 
ing metal abundances very low. Most observed baryons 
must have been inside an early generation of stars, the 
early stars must all have had intermediate initial masses 
in the range 2 — 8Af Q , and they must not have ejected 
much C or N. 

Other astrophysical explanations seem equally unlikely 
(Epstein et al. 1976; Jedamzik & Fuller 1997; Fuller & 
Shi 1997; Famiano, Boyd & Kajino 2001; Pruet, Guiles 
& Fuller 2002; Jedamzik 2002). We can not hide signifi- 
cant D in dust or molecules because neither are abundant 
enough in typical absorption systems. We would require 
> 10% of all the gas to be molecular if HD/H 2 ~ 10" 4 , 
or a large proportion of all heavy elements to be in dust, 
neither of which are likely, except in molecular clouds. 

9.6. Our Estimate for the Primordial D/H from all 

QSOs 

We believe that the best value for primordial D /H is the 
weighted mean (§9.1) of the log D/H values for the first 
five QSOs listed in Table 5: 

logD/H = -4.556 ±0.064, (5) 
(lcr error of 15%), which is equivalent to D/H= 2.78±%?£x 
10 -5 , where the errors are the la errors on the mean, given 
by the standard deviation of the five log D/H values di- 
vided by s/5. We use this error on the mean instead of the 
usual error on the weighted mean because the individual 
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D/H values show more dispersion than we expect. The 
error on the mean depends on just the dispersion of the 
D/H values, and the number of measurements, but not on 
the errors on the values. Had the D/H values been con- 
sistent with a single value (§9.1) we would have used the 
error on the weighted mean, and if the errors on the indi- 
vidual values were also unchanged, this error would have 
been 0.023, or 1/3 of the error we quote. 

The new D/H is 0.6a lower than the value we gave in 
O'Meara et al. (2001), log D/H = -4.52±0.06, because the 
two new values since that paper are both lower. However, 
the error, which is the error on the mean in both cases, 
has not changed significantly. Although the dispersion in 
the D /H values is now larger, the error on the mean is not 
larger because we have two more measurements, Q2206- 
199 and Q1243+3047. 

When we take the primordial D/H from the weighted 
mean, we assumed that the quoted errors on each D/H 
value are too small by the same factor. However, if 
this assumption is not true, other ways of combining 
the measurements of D/H will give a better estimate of 
the primordial D/H. As an example, we could speculate 
that Q1243+3047 and HS 0105+1619 give the two most 
reliable D/H values because Q1243+3047 has the most 
thorough treatment of the errors, and HS 0105+1619 is 
the simplest measurement with the most supporting evi- 
dence. A similar line of argument was explored by Pet- 
tini & Bowen (2001) when they derived a D/H value us- 
ing three DLA systems alone, although the value they 
gave is no longer acceptable because the D/H to Q0347- 
3819 has since increased and now seems the least se- 
cure (§15.2). The best estimate of the primordial D/H 
might now be the weighted mean of the D /H values from 
HS 0105+1619 and Q1243+3047 alone. We might then 
reject the other three measurements because they have 
less data and they were less thoroughly evaluated (§11.1). 
This alternative, two-QSO, mean D/H is not very differ- 
ent: logD/H = -4.604+ 0.032 (7.6% error), or'D/H = 
2.49x 10~ 5 , and the corresponding parameters from SBBN 
are 77= 6.30 + 0.30 x 10~ 10 and Sl b h 2 = 0.0230 + 0.0011, all 
0.8<7 different from the values that we quote in §10, in units 
of the errors in that section. 

10. BBN RELATED COSMOLOGICAL PARAMETERS 

We use SBBN calculations to obtain the r\ and U^h 2 val- 
ues that correspond to the primordial D/H. We use this 
r] value to predict the abundance of the other light nuclei, 
and we compare with measurements. There are differences 
that may be caused by systematic errors. We also compare 
with other estimates of the Qt,h 2 and find good agreement. 

Using the SBBN calculations of Buries, Nollett & Turner 
(2001), our best estimate for primordial D/H leads to 
the following predictions: 7/ = 5.9 ± 0.5 x 10~ 10 , Qf,h 2 = 
0.0214+ 0.0020 (9.3%), Y p = 0.2476+ 0.0010 (predicted 
mass fraction of 4 He), 3 He/H = 1.04 ± 0.06 x 10" 5 and 
7 Li/H = 4.5±o s x 10~ 10 . In the above, the error on Y p 
is the quadratic sum of 0.0009 from the error on the D/H 
measurement and 0.0004 from the uncertainty in the Y p 
for a given 77 (Lopez & Turner 1999). We obtain slightly 
different central values if we use values from Esposito et 
al. (2000a, b). The differences are 10% or less of the error 
from D/H alone, except for 7 Li/H (Esposito et al. equa- 



tions give 4.05 x 10" 10 ) and 3 Hc/H (1.06xl0~ 5 ). In Fig. 
32 we compare the predicted abundances with some recent 
measurements. The vertical band shows the range of 77 and 
rith 2 values that SBBN specifies for our primordial D/H 
value. Measurements of primordial 3 He are consistent, but 
all 7 Li and most 4 He measurements prefer lower 77. 

10.1. Comparison with the Abundances of other light 

nuclei 

Here we discuss the other light nuclei produced during 
BBN - 3 He, 4 He, and 7 Li- and why we believe that D/H 
is preferred over these elements to determine Ot, . 

For several years it appeared that the primordial 3 He 
had not been measured, because of chemical evolution. 
Now Bania et al. (2002) report a limit on the primor- 
dial 3 He/H ratio from their detailed long term study of 
60 Galactic H II regions and 6 planetary nebulae. They 
argue that low mass stars have neither destroyed nor re- 
leased significant 3 He, because 3 He/H changes little with 
metal abundance. In 17 Galactic H II regions for which 
the ionization corrections were relatively simple, they find 
a mean 3 He/H = 1.9 ± 0.6xl0 -5 , which will be an upper 
limit on the primordial ratio if stars have not on aver- 
age destroyed 3 He. They propose that the best value for 
the upper limit on the primordial 3 Hc/H is the value they 
measured for one H II region that has the lowest metal 
abundance in their sample, the third lowest 3 He/H ratio, 
excellent data and a small ionization correction of 22%. 
They then quote 3 Hc/H < 1.1 ± 0.2xl0~ 5 that is consis- 
tent with the value predicted by D/H and SBBN. Given 
the potential complexity of the chemical evolution of 3 Hc, 
the relatively small range and high mean metal abundance 
in the gas where they have made measurements (0.1 < 
[O/H] < —0.5), and the other possible ways of extracting 
the primordial abundance from the data, we suspect that 
the errors are larger than quoted, as with D and 4 Hc. 

The main isotope of He, 4 Hc, is measured in many tens 
of cxtragalactic H II regions to much higher accuracy than 
cither D or 3 He. But 4 He is fairly insensitive to 77 and the 
differences between the measurements allow a large range 
of 77, probably including the value indicated by D/H. It is 
likely (Skillman et al. 1994; Olive & Skillman 2001) that 
the systematic errors were underestimated for many mea- 
surements. Pagel (2000) stated "systematic errors up to 
about 0.005 are still not excluded" , and "Y p is very proba- 
bly between 0.24 and 0.25" , while Fields et al. (2001) state: 
"the systematic uncertainties in the 4 He abundance make 
it difficult to exclude Y p as high as 0.25" . The upper end of 
this range includes the Y p expected from D/H and SBBN. 
However, Peimbert et al. (2002) argue that the systematic 
errors can be ten times smaller. Compared to other mea- 
surements, Izotov & Thuan (1998) report relatively high 
values for 4 He/H that are closest to the D/H prediction, 
but at Y p = 0.2443 ± 0.0013 they are still approximately 
2cr below our prediction from D/H (in units of their er- 
ror). Thuan & Izotov (2002) estimate that their 1998 Y p 
value may be to small by 0.0005 - 0.0010, which still leaves 
Y p 1.5cr below the prediction from D/H. Ballantyne, Fer- 
land & Martin (2000) reject nebulae from the Izotov & 
Thuan (1998) sample that may have significant ioniza- 
tion corrections and find a very high Y p = 0.2489 + 0.0030, 
which is consistent with D/H, but they use an unaccept- 
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able decrease in 4 He/H as 0/H increases. Although Izotov 
& Thuan found absorption lines that might explain why 
Olive, Steigman & Skillman (1997) found less 4 He in some 
objects, it is uncertain whether the Y p is as high as indi- 
cated by D/H. 

7 Li also prefers a lower rj than the value given by SBBN 
and D/H. There are many tens of high accuracy measure- 
ments of 7 Li in the atmospheres of metal poor halo stars. 
There is very little scatter, and hence the abundance ratio 
in these stars is well determined at 1 - 2 xlO~ 10 (Ryan et 
al. 2000; Bonifacio & Molaro 1997; Thorburn 1994). With 
the recent measurement of D/H, 3 He and the Vi b h 2 from 
the CMB, it is looking increasingly likely that the 7 Li in 
these stars is not the primordial value, but rather a fac- 
tor of 3 - 4 less (Tytler, Fan & Buries 1996; Charbonnel 
2002). Authors differ on whether this much depletion is 
reasonable. Pinsonncault et al. (2002) quote a primordial 
7 Li/H = 2.5llg gg x 10~ 10 from models of stellar rotational 
mixing and the measurements of Ryan et al. (1999). This 
value is lcr below the D/H prediction. However, Ryan 
et al. (2000) claim that the depletion correction is only 
0.021q q2 m the log, and that D/H is not consistent with 
7 Li/H. In Fig. 32 we show their primordial 7 Li/H value 
that includes all their corrections, including their depletion 
correction. Vangioni-Flam, Coc & Casse (2000) also be- 
lieve that the 7 Li depletion is small, and hence that 7 Li/H 
rather than D/H, gives the best rj. 

We prefer the rj and fl b h 2 values from D/H rather than 
the lower values from 4 He and 7 Li for several reasons. The 
deuterium abundance is more sensitive to i] that the abun- 
dances of the other light nuclei (Fig. 32). The D/H that 
we see is likely the primordial value (§9) and does not 
need corrections D creation or destruction. No ionization 
corrections are needed. Over the last eight years we have 
established that deuterium is seen in the spectra of a few 
QSOs, and D/H can be measured to approximately 10% 
accuracy in favorable cases. The excess dispersion in the 
D /H values can be explained if some of the measurement 
errors have been underestimated. We account for this dis- 
persion in our estimates of the error on the primordial D /H 
value, the rj value and the other deduced parameters. The 
dispersion is not enough to include the lower rj values pre- 
ferred by 4 He and 7 Li. The D/H value from Q1243+3047 
is similar to that from HS 0105+1619 and the mean D/H 
towards the QSOs is very similar to the primordial D/H 
we expect using a simple model of with Galactic chemical 
evolution to estimate the depletion in the LISM (§9.4). 

Our confidence in the rj and VL b h 2 from D/H is increased 
by the agreement with the fl b h 2 from the anisotropy of the 
CMB, but this does not lessen the need for improved pri- 
mordial abundance measurements. Systematic errors are 
important, or dominant, in most measurement methods, 
and many error terms are poorly known because they are 
hard to estimate. 

10.2. Comparison with other measurements of the baryon 

density 

Aside from using the light element abundances, the cos- 
mological baryon density has been measured in a variety 
of other ways. These methods include the amount of H 
in the intergalactic medium at redshift z ~ 3, the fraction 
of baryons in clusters of galaxies, and most recently and 



with the most accuracy, the variation of the anisotropy of 
the cosmic microwave background on angles of under one 
degree. 

The estimates of ft b h 2 from different CMB experiments, 
listed in Table 7, are consistent with each other and with 
the n b h 2 from D/H and SBBN. The first results from 
BOOMERANG (de Bernardis et al. 2000) indicated a 
much larger value of fl b h 2 ~ 0.036 ± 0.005, however, Net- 
terfield et al. (2002) have since used revised the point- 
ing solutions and include data on smaller angular scales. 
They find n b h 2 = 0.022l°;°° 4 using "weak priors" which 
constrain the age of the universe to > 10 Gyr and the 
Hubble constant to 0.45 < h < 0.9. For the same instru- 
ment, de Bernardis et al. (2002) find Q, b h 2 = 0.022±°;°° 4 , 

0.020lo;o°4 and 0.019tg;°°4 in three complementary calcu- 
lations. Analyzing the measurements from the first season 
of observations with the Degree Angular Scale Interfer- 
ometer (DASI), Prykc et al. (2002) report a very similar 
value, with similar precision: fl b h 2 = 0.0221q'qq3. The 
results from a third experiment, MAXIM A- 1 (Lee et al. 
2001; Stompor et al. 2001) are consistent. Early measure- 
ment on the CMB anisotropy in smaller angular scales by 
the Cosmic Background Imager (CBI) had a maximum 
likelihood with a lower il b h 2 — 0.009, and this likelihood 
dropped by a factor of two for Q, b h 2 = 0.019, and a fac- 
tor of 3 for fl b h 2 = 0.03. Later measurements covering 
40 square degrees on the sky give much more precision. 
Sievers et al. (2002) include the COBE-DMR results and 

find 9, b h 2 = 0.022t£oo9- 

Other methods of measuring the tt b have lower accu- 
racy. The Lyct forest at redshifts z ~ 3 typically indicates 
higher fl b h 2 values (e.g. Hue et al. 2002), while the baryon 
fraction in clusters of galaxies gives consistent values. For 
example, Steigman (2002) multiplied the fl m derived from 
Snla, assuming the universe is flat, by the baryon fraction 
in clusters of galaxies to obtain tt b h 2 = O.OWl^QQg. 

The relevance of D/H measurements is changing. Today 
SBBN and D/H gives the best estimate of tt b h 2 . However, 
as CMB measurements improve, the Q b h 2 from the CMB 
will be as accurate as that from the SBBN. When we use 
n b h 2 from the CMB in SBBN, the D/H is predicted with 
no free parameters, and hence the main value of D/H will 
become to test the physics in SBBN (Kaplinghat & Turner 
2001; Cyburt, Fields & Olive 2001; Steigman, Kneller & 
Zcntncr 2002; Abazajian 2002). Such tests can be made 
now, comparing the abundances of the light nuclei, but the 
measurement errors are not well established, and hence the 
precision will improve when Q, b h 2 comes from the CMB 
and we can use D/H alone to test the physics in SBBN. 

11. SUMMARY 

In this paper, we have presented the detection of D 
towards Q1243+3047. We measured D/H and obtained 
an accuracy nearly as good as the best previous measure- 
ments. The value is slightly lower than the previous mean 
and hence our best estimate for the cosmological baryon 
density from SBBN is slightly higher. 

The most conspicuous absorption system in the optical 
spectrum of Q 1243+3047 makes a strong Lyman break, 
and we found D in the main component of that sys- 
tem. The absorption system has a total column density 
log Nhi= 19.73 + 0.04 cm' 2 in two main components, sep- 
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aratcd by approximately 13 km s _1 . The separation is not 
well determined, because the second component contains 
only approximately 12% of the H I, D I and O I, and hence 
in no case does it make a distinct line. It makes the O I 
slightly asymmetric, and it accounts for a portion of the 
absorption at velocities between the main H and the main 
D, but its only effect on the H I lines is to make them 
slightly wider on the red side. In total the absorption sys- 
tem has at least seven components, but only components 
1 and 2 show O I, and we have seen that nearly all of the 
H I is confined to —40 < v < 40 km s _1 . The other com- 
ponents have log Nhi— 16 cm -2 and with one exception, 
they have no effect on the D/H measurement. 

The column density of the D I in the main component, 
D-l, is well determined, since we see five transitions from 
this absorber, and several of them are well separated from 
other lines, giving accurate 6-values, redshift and log Ndi. 
The column density in the second component, D-2, is less 
well determined because its velocity and 6-value are not 
well known. We explored a variety of models for the rel- 
evant portions of the spectrum, and found the range of 
parameters that gave the lowest x 2 values. Some accept- 
able models have 7% of their Ndi in the D-2 component, 
while others have 28%, the maximum allowed by the spec- 
trum. The range for the column density of D-2 gives the 
range in the total D column density. This range arises 
because the absorption near D-2 can be either D, or H I 
from component H-3 that has log Nhi= 15.90±0.03 cm~ 2 , 
b = 17.0 ± 1.0 km s" 1 and v = -44.3 ± 7 km s" 1 . The 
b- value and v are not very well known for H-3. As the b- 
value increases, or the v decreases, H-3 absorbs more near 
D-2, and the column in D-2 decreases. 

The line near D-l is clearly H or D because it makes a 
Lyman series, and with b = 9.2 ± 0.2 km s _1 , it is much 
narrower than typical H lines. Both the velocity of the 
D-l lines and their width indicate that the line is D and 
not H. Comparison with the width of the O I line gives T 
= 5500 ± 400 K which is cooler than we expect, but there 
is more uncertainty here than indicated by the quoted er- 
ror because the O I and D I are slightly separated in ve- 
locity. We are not as certain as we were for PKS 1937- 
1009, Q1009+2956, Q0130-4021, and HS 0105+1619 that 
we have detected D, because for Q 1243+3047 all the H I 
lines are saturated, and there is extra absorption in other 
components on cither side of the main H I, and hence we 
know little about the velocity distribution of the H I. We 
can not make a detailed comparison the H I and D I ve- 
locity structures that would prove that all of D-l is D. 

While the main error in the D column density comes 
from the uncertain velocity structure, the main error on 
the H column density comes from the uncertain continuum 
shape. The Lya line alone gives the best constraints on 
the log Nhi and it is in the short wavelength side of the 
Lya emission line, where the continuum level is especially 
hard to reconstruct. Had the Lya line been well removed 
from the emission line, the error on the log Nhi might have 
been much less. 

We determined the log Nhi by fitting the continuum, 
Lya forest and the main Lya line simultaneously. We 
model the continuum with smooth curves represented by 
B-splines, and we use an optimizing code to vary the hun- 
dreds of parameters in the model. We re-started the opti- 



mization thousands of times, using different initial values 
for the parameters. All models with the lowest x 2 val- 
ues have log Nhi= 19.73 cm~ 2 , which we use as the best 
fit value. The error on log Nhi is harder to estimate. The 
change in \ 2 values suggests that the la error is near 0.005, 
but we reject this error as too small, because we are un- 
sure if we have adequately explored the parameter space 
or range of models that might be consistent with the spec- 
trum. Instead, we find the range of log Nhi values that 
give models with acceptable \ 2 values, and we represent 
the log Nhi error with a normal distribution chosen to give 
a 20% chance that the log Nhi value is outside the range 
of acceptable models. The errors on the H I and D I are 
then comparable. 

There are now measurements or limits on D/H towards 
seven QSOs, all listed in Table 5. In the appendix, we 
explain why measurements towards two other QSOs are 
no longer considered useful, and why we will not use mea- 
surements towards Q0347-3819, because the D line is not 
seen, and the velocity structure is too uncertain. Ignoring 
the one consistent limit, there are then measurements to 
five QSOs, four from our group. 

These five measurements have a larger dispersion than 
we expect given their quoted measurement errors. There is 
only a 1.5% chance of a larger % 2 value by chance. We sus- 
pect that the dispersion is large because one or more of the 
D/H values is inaccurate, or has errors that have been un- 
derestimated. The work we present here on Q1243+3047 
has reinforced this belief, because we have found that it is 
surprisingly hard to obtain a reliable D/H measurement. 
Many details are relevant (§9.2) and can change the value 
or the error estimate, and we must attempt to model all 
D/H values that might be consistent with the spectrum, 
and not stop with the simplest models. 

The five D/H measurements do not correlate with metal 
abundance, but instead there is a correlation with log Nhi 
that is not expected from cosmology or astrophysics. We 
suspect that this correlation is an artifact, a chance or- 
dering of the D/H values, that themselves show excess 
dispersion. Even though we suspect that some of the D /H 
measurements have larger errors than have been quoted, 
if the relative size of the errors are indicative of the real 
errors, we can use the weighted mean D/H from the five 
QSOs as the best cosmological D/H. For the error, we 
use the error on the mean, and not the weighted mean, 
because the former better represents the dispersion. The 
best cosmological D/H is 0.6cr lower than the value that 
we quoted in O'Mcara et al. (2001) because the two new 
measurements are both below that value. The error is 
unchanged, because the larger number of measurements 
compensates for the increased dispersion. 

The lower value for the cosmological D /H increases the 
tension with the primordial abundances of 4 He and 7 Li, 
both of which prefer lower values for r\ and fl^h 2 . How- 
ever, considering the difficulties in obtaining these primor- 
dial abundances, we consider that the relative abundances 
of the elements are consistent with SBBN. 

For several decades now most measurements have given 
less 4 Hc than we predict, although the dispersion of those 
measurements suggests that this difference is because of 
systematic errors. A recent estimate of the primordial 3 He 
abundance by Bania et al. (2002) agrees exactly with D/H, 
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but they rely on one measurement and we still do not un- 
derstand the change in 3 He caused by Galactic chemical 
evolution. The abundance of 7 Li in halo stars is a fac- 
tor of 3 - 4 less than predicted. Significant 7 Li may have 
been destroyed in these stars, although the amount of de- 
struction implied by SBBN and D /H is near the maximum 
allowed. There is a clear need for more work on all four 
nuclei. 

We believe that Deuterium is the best baryometer 
(Schramm and Turner 1998) among the light nuclei for 
the reason given at the end of §10.1. The agreement be- 
tween the ftth 2 from D/H and that from the CMB adds 
confidence. These two fi^/i 2 measurements use very differ- 
ent physics, at very different epochs in the universe. The 
agreement implies that we understand the main physical 
processes that effect the observations of both SBBN and 
the CMB anisotropies. 

11.1. New methods used in this paper 

For Q 1243+3047, we have attempted to improve the 
methods that we use to measure D/H. We employed a 
more accurate relative flux calibration, a more thorough 
explorations of the models that can represent the D and 
H absorption, more realistic continuum fitting and better 
representation of the Lya forest and other blended ab- 
sorption. The D/H value and the errors should be more 
reliable. 

We have paid much more attention to the accuracy 
of the relative flux calibration of the HIRES spectra of 
Q1243+3047 than we did for any of our past publications. 
We have several low resolution spectra from the Lick Kast 
spectrograph, and also a spectrum from ESI. These spec- 
tra have allowed us to compare several methods of flux 
calibrating our HIRES spectra, and to verify that we get 
similar results if we start with different calibration spectra. 

Our continuum fitting procedure is more general than 
the methods we used for our measurement of D /H in other 
QSOs. For PKS 1937-1009 and Q1009+2956, we fit the 
continuum twice; once as a smooth curve passing through 
the highest flux levels, and again as one smooth curve for 
each portion of the spectrum which was adjusted to give a 
good fit. For HS 0105+1619, the continuum was a single 
smooth curve which passed through the high points in the 
flux which appeared to be least absorbed by the Lya for- 
est. We estimated the errors by moving this curve up and 
down by a factor, such that it was below the flux in many 
places, or above the flux in all places. The B-splinc con- 
tinua that we fit to Q 1243+3047 should give make our Nhi 
estimate more reliable. The B-spline continuum model al- 
lows for a wide variety of smooth shapes, and we are able to 
fit the continuum parameters simultaneously with the ab- 
sorption line parameters. In addition, the B-spline gives us 
the ability localize continuum degrees of freedom in wave- 
length, which we have had difficulty doing in past work. 
These continuum can be adjusted iteratively to allow for 
the effects of line absorption, and hence we no longer have 
to require that the flux is at the continuum level in some 
pixels. 

For Q1243+3047, we fit the continuum, Lya forest and 
the DLA all simultaneously. For HS 0105+1619, we fitted 
the continuum before and independently of the Lya forest 
and the DLA. When we used methods like those that we 



used for HS 0105+1619 on Q1243+3047, we found a sig- 
nificantly lower Nhi value, probably because that method 
required that the flux in at least a few pixels is at the 
continuum level. We have seen here that the continuum 
level, the Lya forest and the DLA are connected in a com- 
plex way. In contrast to HS 0105+1619, where we may 
have use too few degrees of freedom to model the contin- 
uum, for PKS 1937-1009 and Q1009+2956, we may have 
allowed too many degrees of freedom. 

For Q1243+3047, we have also attempted a more thor- 
ough exploration of the models that might represent the 
spectrum. For PKS 1937-1009 and Q1009+2956, we used 
a minimization method to find the optimum parameters. 
Here, for the log Ndi, we use a full grid search, where each 
parameter takes on every allowed value, ensuring a thor- 
ough search. 

For Q 1243+3047, we have also used different ways of es- 
timating the measurement errors. For PKS 1937-1009 and 
Q1009+2956, we used the 8\ 2 method that we used here 
for the log N D i alone. For HS 0105+1619, we quoted the 
error on log Nhi of 0.009 cm -2 , compared to 0.04 cm -2 
for Q1243+3047. The smaller error for HS 0105+1619 is 
the weighted mean of three estimates of log Nhi from dif- 
ferent parts of the Lya line. It is reasonable that the value 
is smaller than for Q1243+3047 because the DLA suffers 
less contamination near its core. To estimate the error 
from the continuum, we shift the continuum up or down 
by some factor, but we did not change its shape. 

11.2. Five Outstanding Issues 

We end with six outstanding issues. Some are specific to 
our D/H measurement from Q1243+3047 and others are 
long standing. 

Velocity offsets. The O I absorption is centered 
2.8 km s _1 away from the center of the D absorption, not 
enough to change the identification of the D, but enough to 
raise concerns that we may have misunderstood the veloc- 
ity structure of the absorption system (§5.1). The deepest 
components of C II and S II are also offset from O-l by 
2 km s , in the other direction from the D-l component. 
These offsets make it to estimate the column densities of 
different ions in the same gas, and hence the ionization is 
uncertain. 

Missing D components. The main error on the 
log Ndi comes from the amount of D in the component 
near 13 km s _1 . This component is hard to recognize and 
fit because it is weak and blended with the main compo- 
nent (Fig. 6). There might be D in other components 
that we have mistakenly modeled as H (§4.2.6). We have 
not included this possibility in the log Ndi error because 
we do not know how to estimate the chance that it has 
occurred. This additional D could be >10% of the total if 
these components have [O/H] < —3. 

Accuracy of the log Nhi error. We spent consid- 
erable time investigating and debating the error associ- 
ated with log Nhi (§6). We mentioned three such error 
estimates in this paper: the covariance matrix from the 
optimized fits, the Sx 2 method, and the range of log Nhi 
that gave fits with acceptable x 2 ■ They gave very differ- 
ent errors: 0.0002, 0.005 and 0.04 cm -2 respectively. We 
know that our estimate from the covariance matrix is un- 
reliable, and we mention it only because it is the standard 
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error quoted in most QSO absorption line studies. 

The error that we use, from the range of acceptable fits, 
is unconventional. It could be in error, perhaps by a fac- 
tor of two, for several reasons. The largest acceptable x 
value might be incorrect (§6.5.4), we may have been over 
generous in allowing Lya forest lines with b = 150 km s _1 
(§6.5.2), or we may have considered too narrow or too 
wide a range of Lya forest (§6.5.2) or continuum shapes 
(§6.5.3). 

Dispersion in the D/H values. While we suspect 
that the excess dispersion arises from underestimated mea- 
surement errors, and we discuss circumstantial evidence 
for this (§9.2, 9.3, & 11.1) we will not know whether this 
is correct until more reliable measurements are made. 

We are struck that the D/H values and their errors de- 
pend more on the choice of the models that we use to 
estimate D/H, than upon the details of the fit to the data 
and the x 2 value (§9.2). A vital role of high S/N spectra 
is to guide us to examine an appropriate set of models. 
The errors on D/H values are hard to estimate because it 
is hard to know whether we have explored a full range of 
models. 

Correlation of D/H with log Nhi- The correlation 
(Fig. 31) looks rather impressive, because all five mea- 
surements, and the limit, allow a monotonic decrease of 
D/H with increasing log Nhi- We suspect the correlation 
is spurious, and a random accident of the measurement er- 
rors. We do not know of any single systematic error that 
could explain all of the effect (§9.5). Rather, a variety of 
errors may be involved. 

The Lack of Precise Agreement of the Measured 
Abundances in SBBN. The measurements of the light 
nuclei abundances have preferred slightly different r\ val- 
ues since our first measurement of D/H (Tytler, Fan & 
Buries 1996). Although Izotov & Thuan (1998) reported 
higher Y p values that were closer to the predictions from 
SBBN and D/H, our D/H values have decreased over the 
years, exacerbating the difference in the favored rj values. 
There are excellent reasons to believe that measurement 
problems account for the lack of precise agreement. 
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13. APPENDIX I: B-SPLINE CONTINUUM MODEL 

We model the continuum level with a B-spline curve. 
A B-spline is defined by a set of control points, 



Po, Pi, . . . , P n , n > 3 where n + 1 is the number of control 
points defining the curve. The B-splinc passes near to the 
control points, but does not generally pass through them. 
A B-spline is defined by the control points as a series of 
n — 2 piecewise parametric polynomial segments. The z'th 
segment is given by 
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where t is the parametric variable and ranges from to 
1. The full B-spline curve, including the joints between 
segments, has continuity of order two. That is the curve 
itself, and it's first and second derivatives are continuous. 

In the past we have defined continuum levels with 
Chebyshev polynomials and with cubic splines, which are 
the options available in the IRAF task continuum. We de- 
cided against using such representations of the continuum 
in this work because they both define global curves - if 
you change one control point in a cubic spline or a single 
coefficient in a polynomial the entire curve is affected. In 
contrast, when Pi is changed, only the B-spline segments 
between Ci-2 (t) and Cj+i (i) are affected. This locality has 
two advantages. First, it is easy to define the continuum 
over a large wavelength range in a piecewise fashion e.g. 
modifications to the emission line continuum do not mess 
up the continuum level defined near the Lyman limit. Sec- 
ondly, it greatly simplifies the task of efficiently computing 
the continuum changes that occur as a result of changing 
a single control point, as often occurs during non- linear 
optimizations. 

14. APPENDIX II: OPTIMIZER 

Here we describe the code we used to vary the parame- 
ters of a model for the spectrum. In our previous absorp- 
tion line work, we have usually fit absorption line parame- 
ters using the \ 2 optimize VPFIT, which was kindly made 
available to us by R.F. Carswcll and J. Webb and collabo- 
rators. We were unable to use VPFIT for this work because 
we desired to optimize both our continuum and absorption 
line models simultaneously. In addition, we wanted to be 
able to optimize many hundreds of parameters simultane- 
ously, which we had difficulty doing with standard tools. 
The optimizer we developed to address these issues imple- 
ments the Levenberg-Marquardt algorithm where we reg- 
ularize the normal equations by inverting the curvature 
matrix by means of a truncated singular value decompo- 
sition (SVD). 

Our algorithm is based on the Press et al. (1992) rou- 
tine mrqmin, and we refer the reader to Press et al. for a 
general description of the Levenberg-Marquardt (LM) al- 
gorithm, which we will not repeat here. We only describe 
the differences between our algorithm and the algorithm 
implemented by mrqmin. 

Using the Press et al. notation, a key step in the LM al- 
gorithm is to solve the following set of linear equations for 
the changes Sai that need to be applied to the parameter 
vector ai during the next iteration of the algorithm 



M 



^ a k iSai = /3 k 



(7) 



2 = 1 
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where M is size of the parameter vector. The curvature 
matrix a is defined as 

1 d 2 X 2 

a kl = n q_ o_ ( 8 ) 



and /3 is given by 



2 daedal 



0k = --^ 
2 9aj. 



(9) 



In the most implementations of the LM algorithm, Equa- 
tion 7 is solved by direct inversion of the curvature matrix. 
For example, mrqmin does this by Gauss- Jordan elimina- 
tion. 

The problem with a direct solution of Equation 7 is that 
when fitting absorption lines the curvature matrix is fre- 
quently ill-conditioned. When a has high condition num- 
ber, a direct inversion gives a solution vector Sai that is 
usually nonsense. As a result, the LM algorithm will fail 
to find acceptable parameter updates and will stop. 

Unfortunately, it is very easy for to get an ill-conditioned 
a during an absorption line optimization - all that is 
required is a pair of nearly degenerate parameters. In 
practice this seems to occur when fitting several heavily 
blended lines, or when both broad lines and continuum de- 
grees of freedom arc in the parameter set. There are likely 
many other scenarios which result in an ill-conditioned a, 
as we have observed that the LM algorithm will almost 
always have a few iterations with an ill-conditioned a if 
the number of optimized parameters is large enough. The 
critical number of parameters appears to be something in 
the range of 30-50 though we have not investigated this 
carefully. We just note that in the course of this work an 
ill-conditioned a occurred frequently enough to make the 
standard LM algorithm nearly useless. 

In most situations when the curvature matrix is ill- 
conditioned, there are still parameter updates Sai available 
that would significantly improve the \ 2 between the model 
being optimized and the underlying data. To find these, 
and to thus improve the performance of our optimizer, our 
algorithm regularizes ciki before solving Equation 7. We do 
this by inverting ctki via a truncated singular value decom- 
position. Briefly, we compute the SVD of the a using the 
LAPACK routine dgesdd, which computes orthogonal ma- 
trixes U, V and the diagonal matrix X = Diag(eri, . . . , o~ n ). 
U, V and £ are defined such that 

a = UT,V* (10) 

Because U and V are orthogonal and £ is diagonal it is 
trivial to invert a once it's SVD is known 

or 1 = VIT 1 U* (11) 

which works as long as all of the au are non-zero. The SVD 
(Equation 10) always exists, even if a is ill-conditioned or 
even singular. Unfortunately, if a is ill-conditioned, Equa- 
tion 11 does no better than direct methods for computing 
a -1 . However, the magnitude of the o~k tell us which of 
the orthogonal vectors in U arc responsible for the near 
singular behavior. The standard trick is to compute a -1 
via 

aT 1 =VE + U* (12) 



where 



In this scheme h is called the regularization parameter. 
Setting it is something of a black art. If h is set too small, 
the near singular components of U will corrupt a -1 . If 
h is set too large, otherwise good component vectors of 
U will not be updated, and the optimizer will perform 
poorly. We could find no guidance from the literature on 
an appropriate value for h - the general advice given is too 
experiment and see what works for the application. This 
is exactly what we did, and in our algorithm h = 10~ 8 . 
Our entire algorithm is implemented in double precision. 

This is the algorithm that was used for all of the opti- 
mizations performed in this work. This optimizer works 
very well in practice, and we have routinely performed 
successful optimizations of over 500 parameters, many of 
which were highly blended Lya forest lines. 

14.1. Error estimates 

The standard way to estimate the errors in the fitted 
parameters is to estimate the covariancc matrix C from 
the curvature matrix a 



C = a 



(14) 



E+ = Diag(a+), 



+ _ 1 1 / (T k ifofc > h 
k 10 otherwise. 



(13) 



The parameter errors are then just given by the square 
root of the diagonals of C. 

This is exactly the error estimate that we use in our 
optimizer. We directly invert a using Gauss-Jordan elim- 
ination instead of using the truncated SVD regularization 
described above. This is because the validity of an ap- 
proximate inverse is easily tested when updating parame- 
ter values: if the solution to Equation 7 results in a better 
X 2 , the solution has some value. However, we are not sure 
how an approximate inverse affects the errors inferred from 
the covariance matrix, so we just use a direct method. 

Note that the estimated errors are not likely to be valid 
when a is ill-conditioned. In practice, they generally ap- 
pear to be much to small when the curvature matrix is ill- 
conditioned. We did a small amount of Monte-Carlo test- 
ing which indicated that the estimated covariance matrix 
appeared to be accurate when optimizing a small number 
of lines with a well conditioned covariance matrix. How- 
ever, it appeared from the same testing that the estimated 
covariance matrix begins to go bad long before a becomes 
ill-conditioned enough to jeopardize the solution of Equa- 
tion 7. Because we are sure that Equation 14 is a valid 
estimate in only very simple optimizations, we use error 
estimates from Equation 14 only cautiously in this work. 

15. APPENDIX III: DISCUSSION OF OTHER d/h 
MEASUREMENTS 

In this appendix we discuss D /H measurements and lim- 
its to other QSOs. In some cases the absorption systems 
are too complex, or existing spectra are inadequate to give 
convincing measurements or constraints. Related issues 
are discussed in reviews by Tytler & Buries (1997), Tytler 
et al. (2000) and Lemoine et al. (1999). 

Towards Q0014+813 Songaila et al. (1994) reported an 
upper limit of D/H < 25 x 10~ 5 in the z abs = 3.32 LLS. 
Using different spectra, Carswell et al. (1994) reported 
< 60 x 10 -5 in the same object, and they found no rea- 
son to think that the deuterium abundance might be as 
high as their limit. Improved spectra (Buries, Kirkman & 
Tytler 1999) showed that D/H < 35 x 10" 5 . A high D/H 
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value is allowed, but is highly unlikely because the absorp- 
tion near D is at the wrong velocity, by 17 + 2 km s _1 , it is 
too wide, and it does not have the expected distribution of 
absorption in velocity, which is given by the H absorption. 
Instead this absorption is readily explained entirely by H 
(D/H ~ 0) at a different redshift. 

Towards PG 1718+4807 the LLS at z abs = 0.701 was 
also considered as a possible site showing a large D/H value 
(Webb et al. 1997), but others argued that this was not the 
most convincing interpretation of the HST spectra (Lev- 
shakov, Kegel, & Takahara 1998; Tytler et al. 1999). New 
higher resolution spectra from HST (Kirkman et al. 2001) 
show the velocity structure of the H for the first time. We 
found that it is unlikely that this QSO gives any useful 
information on D/H, because the absorption near D is at 
the wrong velocity and has the wrong line width to be all 
D. Hence it is probably mostly contaminating H. 

Molaro et al. (1999) claimed that the z=3.514 absorp- 
tion system towards AMP 08279+5255 showed low D/H, 
but they and Levshakov et al. (2000) note that since only 
the Lya line has been observed, the hydrogen velocity 
structure and the H I column density are poorly known, 
and the D feature can be fit using H alone. 

15.1. Q 2206-0 199 

Pcttini and Bowen (2001, hereafter PB) report D/H 
= 1.65±0.35xl0" 5 in the z = 2.076 DLA towards Q2206- 
199. Since the D is measured in a very low S/N HST spec- 
trum it is hard to determine whether the errors are rea- 
sonable, and whether D is seen. The reported D/H value 
is very low compared to the D/H we expect from Galac- 
tic chemical evolution and the D/H in the LISM (§9.4). 
We now discuss some ways in which the errors might have 
been underestimated. 

We concentrate on the Ndi which seems to be less re- 
liable than the log Nhi- The authors solve for the Ndi 
by fitting features at the position of D in three Lyman 
series transitions: Ly-7, Ly-9, and Ly-12. The Ndi value 
depends on the 6- value chosen for the D. The D lines are 
too noisy to give this 6-value, so the authors interpolate 
between the widths they measure for the 6-values of the 
H I and the metal lines, giving 10.6 km s _1 . We ques- 
tion whether this 6-value is accurate, perhaps because the 
6-value of the H I is itself uncertain. 

Although the low ion transitions such as Fe II and Al II 
indicate a single component nature for the gas, Prochaska 
& Wolfe (1997) also fit multi-component CIV and SilV 
gas on either side of the v = position of the H which 
shows D. These additional components may have enough 
H absorption to widen the Lyman lines that were used to 
get the 6-value of the main H I. As such, we explore the 
implications of a lower 6(D) on the range of D/H allowed 
by the data. If we assume the gas is cold, in the range 
of 300-400 K, and motions are dominated by turbulent 
velocities, then values of log Ndi giving D/H ~ 2.9xl0~ 5 
are consistent with the 1 a errors obtained by PB's fit to 
Ly-9, but not to Ly-12. 

We consider the other sources of error relevant to the 
determination of log Ndi • In the Lyman limit region pre- 
sented by PB, very few pixels reach the continuum in 
their normalized data. In particular, for Ly-9 and Ly-12, 
no pixels reach the continuum within v ~ 150 km s -1 of 



the position of D. Due to the high level of contamination 
throughout the Lyman limit, the continuum is very poorly 
constrained, and could contribute significant errors to the 
log Ndi- 

Continuum placement aside, contamination also affects 
the determination of log Ndi in this system. In particu- 
lar we note that the fit to the contamination on the large 
wavelength side of the D feature in Ly-9 seems incompat- 
ible with the size of the 1 a errors given to the D line. 
Since this additional absorption adds significant optical 
depth at the position of D, it might change the log Ndi by 
more than 2a. 

Our main concern with this D/H measurement is that 
the HST spectra have low S /N and provide very little sup- 
porting evidence. However, the ions do show a simple ve- 
locity structure, and hence the feature may be D with the 
measured D/H. 

15.2. Q0347-3819 

D'Odorico, Dessauges-Zavadsky k Molaro (DDM, 2001) 
report D/H = 2.24 ± 0.67x 10~ 5 in the z=3.025 absorp- 
tion system to Q0347-3819, while Levshakov, Dessauges- 
Zavadsky, D'Odorico & Molaro (LDDM, 2002) report a 
very different value, D/H = 3.75 ± 0.25xl0 -5 , from a so- 
phisticated analysis of the same spectrum. We are not 
convinced that D has been detected in the spectrum of 
this QSO. The D lines are completely blended with the H 
lines and in no case can we see D clearly separated from 
the H. Hence we do not know the central velocity, velocity 
structure or 6-value for the absorption near D. This ab- 
sorption could be partly or predominantly H. It is hard to 
obtain either a secure value or a limit on D/H, because 
the velocity structure for this absorption system is critical 
and not well enough known. We will discuss this system 
in detail here because many of the issues are very familiar 
from our discussion of Q1243+3047. 

The total H I column should be well determined from 
the damping wings of the Lya line. The fits shown in 
Fig. 11 of LDDM seem to over-absorb in Lya near —2800 
km s _1 , and from +450 to +750 km s _1 for Ly/3 which 
suggests that the model has excess log Nhi or that there 
might be errors in the flux calibration or continuum level. 

The metal lines show many components that can con- 
tribute to the Lyman scries H I lines. Several components 
may have enough log Nhi to also show D. The low ions, 
especially H2, indicate the velocity for much of the H I, 
and the expected D. There is absorption on the blue side 
of Lyman lines 8, 10 and 12 at the position expected for 
D, which resembles a Lyman series, and hence it may be 
caused by D I or H I or a combination of the two. 

To determine the total log Ndi that is associated with 
the total log Nhi we need to know the velocity structure of 
the D and the H. DDM assumed that the H I is at the ve- 
locities of three of the components seen in the metal lines, 
and they assumed that the H I columns are distributed 
approximately like the metals. They fit the H I lines by 
adjusting the log Nhi and 6-values of these components. 
LDDM stated that their model may not be unique, which 
seems likely because their model seems to have more de- 
grees of freedom than are needed to fit the data. We expect 
that two components might suffice, one to fit each side of 
the H Lyman series lines. 
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It is difficult to measure the amount of H I in the differ- 
ent components because their lines are all blended. This is 
especially true for the main component (v = 0, component 
3 in DDM) because in no case do we see H I absorption 
from this component alone. The other components, which 
are on either side of it, account for both sides of the higher 
order Lyman lines. 

DDM find that additional absorption is needed near the 
velocity expected for D. It is difficult to measure a col- 
umn density for this extra absorption, whether of not it is 
D, because we do not know its velocity structure. LDDM 
demonstrate that the fit tabulated by DDM, and hence 
the D/H value, is not unique, and that the differences arc 
highly significant because some of the H I, and we expect 
a similar proportion of the D I, is in a saturated compo- 
nent with an unusually low b ~ 2.9 km s _1 . The log Ndi 
is sensitive to the b- value of the D, and the proportion of 
the H I and D I that is in this component, neither of which 
can be measured. 

The log Ndi also depends on the amount of H I absorp- 
tion on the red of the D position. The more flux absorbed 
by the H I, the less can be absorbed by the D, and hence 
the lower the D/H. The total log Nhi is not changed be- 
cause it is fixed by the damped Lya line that is insensitive 
to the velocity structure. If we adjust the relative amount 
of H I in the two main components, both the log Nhi and 
the 6-values, which are not known, then the amount of 
absorption that could be D also changes. There is also 
a component near —65 km s _1 , seen in C IV and Si IV 
(Prochaska & Wolfe 1999), which might also have enough 
H I to effect the absorption on the red side of where D 
is expected in Ly-8, 10 and 12. The uncertainty in con- 
tinuum near Ly-10 and 12, two of the three regions used 
to get log Ndi, and the lines blended on the blue side of 
the expected positions of the D lines further increase the 
uncertainty in the log Ndi- 

DDM and LDDM get very different log Ndi and hence 
D/H because they make different assumptions about the 
velocity structure. They both assume that the H I, and 
the D I, have the components at the velocities that show 
low ionization metal lines, and they assume that the pro- 



portion of the H I in these components is similar to that 
of the metals, which will not work if the metal abundances 
vary. DDM find that the D can be fit with two main com- 
ponents with b- values of 14.1 ± 0.5 and 16.2 ± 3.0 km s _1 . 
LDDM assume that most of the D has b ~ 2.8 km s _1 
which produces saturated D lines that require much larger 
log Ndi to explain a given amount of absorption. It is rea- 
sonable that much of the H I and hence the expected D, 
has b closer to 2.8 than to 15 km s _1 , and hence the D/H 
will be larger than reported by DDM. However we do not 
know how much of the H I has this low b- value, and hence 
we do not know whether the D/H is lower, similar to, or 
high than proposed by LDDM. 

LDDM tie the unknown velocity structure of the H I 
and D I to that of the H2 and various metal lines. This is 
not a unique procedure because in this system the metal 
lines themselves differ: Fe II (1143.2, 1133.7, 1125.4) has 
different velocity structure from Si II 1808 and other ions. 
LDDM note that the Fe II lines may lack column in the 
main component at v = because the Fe there is locked 
up in dust. 

The metal abundance in the gas where D would be seen 
is uncertain because this gas has a very low 6-valuc and 
there may be dust. LDDM find [Zn/H] = -1.51 ± 0.05, 
but elements C, P and Ar are higher by 0.3 dcx, while 
O and Si are higher by 0.7 dex. These enhancements are 
unexpected because Zn should be the least depleted onto 
dust, and we do not expect such a large enhancement of O 
and Si. Although LDDM quote [Si/H] = -0.95 ±0.02, we 
expect that the uncertainty is much larger than 0.02 dex, 
and we anticipate that further analysis will show that the 
total abundances including metals in dust, is lower than 
this. 

We prefer not to take the DDM D/H value as a lower 
limit, because we do not know whether D has been seen. 
Before we could take the LDDM value as an upper limit 
we would need to know that there is at least as much H 
absorption near the D as they assume, and that the pro- 
portion of the D in the low b- value component is at least 
as much as they assume. This absorption system is not 
ideal for the measurement of D/H. 
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Table 1 

OBSERVATIONS OF Q1243+3047 



Instrument 


L/dLc 


t 4-" re- 
integration I mic 




Wavelengths 


Observed 






(seconds) 


(arcsec) 


(A) 


Kast 


February 13, 1997 


3300 


3 


3120 - 


5800 


Kast 


May 10, 1999 


6982 


2 


3137 - 


7188 


Kast 


May 11, 1999 


7200 


2 


3137 - 


7188 


Kast 


May 17, 2001 


5400 


2 


3191 - 


5881 


Kast 


May 17, 2001 


5400 


2 


3191 - 


5881 


ESI 


January 11, 2000 


1364 


1 


4000 - 10,000 


HIRES 


April 15, 1999 


3600 


1.14 


3494 - 


5842 


HIRES 


April 16, 1999 


7200 


1.14 


3168 - 


4705 


HIRES 


April 17, 1999 


8100 


1.14 


3168- 


4705 


HIRES 


April 17, 1999 


8100 


1.14 


3168- 


4705 


HIRES 


March 12, 2000 


9000 


1.14 


3214 - 


4705 


HIRES 


March 12, 2000 


9000 


1.14 


3214 - 


4705 


HIRES 


March 13, 2000 


7200 


1.14 


3214 - 


4705 


HIRES 


March 13, 2000 


7200 


1.14 


3214 - 


4705 



a We list the local calendar date at sunset at the start of the night. 



Table 2 

RESOLUTION AND S/N OF SPECTRA 



Spectrograph 


Slit 


Pixel Width 


S/N a 


S/N a 


FWHM 




(arcsec) 


(kms -1 ) 


(3250A) 


(4250A) 


(kins -1 ) 


Kast 


2 


105 b 


20 


60 


283 ± 25 


ESI 


1 


11.5 




40 


63.2 ±3.0 


HIRES 


1.14 


2.1 


10 


90 


8.0 ±0.2 



a S/N per pixel. 

''Mean value. We measure variation with wavelength and from spectrum to spectrum. 



Table 3 

IONS IN THE z ~ 2.526 ABSORPTION SYSTEM TOWARDS Q1243+3047 



Ion 


log N 
(cm ) 


b 

firm C~ 1% 1 
I Kill a J 




z 


V 

I Kill b 1 


H I total* 


19.73 ±0.04 










H I (H-3) 


15.90 ±0.03 


17.0 ± 1.0 


2 


.525171 


-44.3 ±7.0 


H I (H-l) 


19.63 


14.8 ±2.9 


2 


,525659 


-2.8 


H I (H-2) 


19.05 


10.9 ±3.3 


2 


,525804 


9.5 


H I (H-4) 


16.25 ±0.02 


25.8 ±0.9 


2 


,526939 


106.0 ±0.7 


H I (H-5) 


16.35 ±0.02 


26.6 ±0.5 


2 


,528108 


205.4 ±0.3 


D I total 


15.113 q 










D I (D-l) 


15.058 ±0.03 


9.2 ±0.2 


2 


,525659 


-2.8 ±0.6 


D I (D-2) 


14.191 ±0.10 


< 12 


2 


,525804 


9.2+? «- 

— 3.5 


I (0-l) b 


13.570 


6.766 


2 


525692 


o.o c 


O T (D-9\ h 


1 9 755 


6 766 

U . i UU 


2 


525848 


1 3 3 C 


Si TT 

Ol 11 


1 1 884 

1 1 .004: 


3 44 


9 


^9^900 


—41 


Si TT 

Ol 11 


1 1 783 

11.1 oo 


O.Oi/ 




^9^349 


— 9Q 8 


Si TT 

Ol 11 


1 9 090 


7 8^ 
i .00 


9 


O^OOOO 


lu.U 


Si TT 

Ol 11 


1 1 830 

1 1 .OlJU 


3 7fi 


2 


O AOxjOO 




Si TT 

Ol 11 


1 9 766 


4.46 


2 


OAO t z*o 


3 1 c 


Si TT 

Ol 11 


1 9 473 


5 QQ 

tj . act 


9 




1 4 C 

1U.4: 


Si TT 

Ol 11 


1 1 ^08 

1 1 . ooo 


97 


9 




94 & 


Si TT 

Ol 11 


1 9 1 90 

1 £j . 1 ZiU 


7 41 

1 .4:1 


9 


^9fi0&3 


33 3 


Si TT 

Ol 11 


1 1 7M 
ii.iji 


5 Q7 


9 




QQ 


C II 


12.521 


3.03 


2 


52521 5 


—40.5 


C II 


12.987 


11.06 


2 


525324 


—31.2 


C II 


12.693 


9.97 


2 


,525552 


-11.9 


C II 


12.956 


7.58 


2 


,525626 


-5.6 


C II 


13.344 


5.21 


2 


,525708 


1.4 C 


C II 


13.260 


7.17 


2 


,525872 


15.4 C 


C II 


12.921 


12.73 


2 


,525959 


22.8 


C II 


13.150 


10.40 


2 


,526076 


32.7 


C II 


12.916 


10.59 


2 


,526834 


97.1 


Al II d 


< 12.2 


7 






C 


Si III C 


< 13.0 


7 






C 


C IV f 


< 12.75 


7 






C 


C IV 


13.190 


38.81 


2 


,525750 


4.9 


C IV 


13.256 


20.98 


2 


,526171 


40.8 


C IV 


13.248 


29.53 


2 


,526985 


110.0 


C IV 


12.514 


14.11 


2 


,528017 


197.7 


Si IV f 


< 12.4 


7 






C 


Si IV 


12.905 


37.64 


2 


,525560 


-11.2 


Si IV 


12.730 


17.21 


2 


,526145 


38.6 


Si IV 


12.709 


24.33 


2 


,526944 


106.5 


Si IV 


11.920 


10.71 


2 


,527978 


194.4 



a Includcs components H-l and H-2 only. 

^Values correspond to model where the 6 of O-l and 0-2 have been set to 
equal each other. 

c The column density for this component was used to constrain the ionization 
in components 1 & 2. 

^This ion was observed in the ESI spectrum only. The quoted upper limit 
on the column density is obtained when we fix b = 7.0 km s _1 , v = km s _1 
and fit the entire region from -60 km s _1 to +60 km s — . 

e Not fit because strongly blended. We quote the maximum column near 
v = for 6 = 7 km s that is consistent with the residual flux. 

The upper limit comes from a fit with fixed 6 = 7 km s~ 1 and v = km s" 1 . 
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Table 4 

INFERRED PHYSICAL CONDITIONS WHERE D/H IS MEASURED 



QSO 


logtf I/H 


log J? 


Size 


Hydrogen Gas Mass 






(cm~ 2 ) 


(kpc) 


(M ) 


PKS 1937~1009 a ' b 


-2.35,-2.29 


20.05,19.74 


0.9, 0.4 


3.9 x 10 5 ,2.9 x 10 4 


Q1009+299 b < c 


-2.97,-2.84 


19.90, 19.93 


1.8, 1.5 


1.1 x 10 6 ,7.5 x 10 6 


HS 0105+1619 b < d 


> -0.1 


< 19.52 


< 1.1 


< 1.6 x 10 5 


Q1243+3047 6 


-0.69 


20.42 


2.7 


8.0 x 10 6 


Q2206-199 f 


unknown 6 


unknown 


unknown 


unknown 


Q0347-3819 h 


> -0.3 


< 20.93 


0.014 


211 


0:0130-403' 


-3.4 


20.06 


30 j 


4.5 x 10 8 



a We list the parameters for each of the two components, where available, from Tytler, Fan & 
Buries (1996); Buries & Tytler (1998a). 

^Scaled from previous references to J912 = 10 -21 ergs cm -2 s — 1 Hz" 1 sr — 1 . 

c We list the parameters for each of the two components, where available, from Tytler & Buries 
(1997) and Buries & Tytler (1998b). 

d O'Meara et al. 2001. 
c This paper. 

f Pettini & Bowen 2001. 

s Currently, there is not enough diagnostic information to measure the ionization state in this 
absorber, yet it is presumably neutral, given the log Nhi- 

h Levshakov et al. 2002. 

'Kirkman et al. 1999. 

J Calculated from log nu = —2.9 cm -3 from Kirkman et al. 2000 and log nui/nu = —3.4 from 
O'Mcara et al. 2001 . 
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Table 5 

D/H MEASUREMENTS TOWARDS QSOS 



QSO 


ZDH 


D/H ±lcr 


logD/H 


A, a 6(D) (km s- 1 ) 






(xlO~ 5 ) 




predicted observed 



PKS 1937-1009 b 


3.572 


3.25 ±0.3 


-4.49 ±0.04 


±1.65 


12.5 ± 2.1 c 


14.0 ± 1.0 


Q1009+299 d 


2.504 


3.98«;f 7 


-4 4n+ 006 


±1.95 


13.5±0.5 C 


15.7 ±2.1 


HS 0105+1619° 


2.536 


2.54 ±0.23 


-4.596 ± 0.040 


-1.00 


10.1 ±0.3 C 


9.85 ±0.42 


Q1243+3047 f 


2.525675 


/1O+0.35 


-4 61 7+ 058 

1-OJ-' -0.048 


-1.05 


11.3 ± 1.8 


9.2 ±0.2 


Q2206-199 g 


2.0762 


1.65±0.35 h 


x yo+0.08 


-2.17 


10.6 




Q0347-3819 1 


3.024855 


3.75 ± 0.25 h 


-4.43 ± 0.03 


±3.35 


3,14.1,16.2 




Q0130-403 j 


2.799 


< 6.8 


< -4.17 




16.2 ±0.3 C 





a Xi = (Yi — mean)/cr(Yi), where Yj = log(D/ H)i and we use the weighted mean of the first five QSOs, 
log D/H = -4.556 ± 0.064. 

b We list combined results for the two components, from Tytler, Fan it Buries (1996); Buries & Tytler 
(1998a). 

c Calculated from the published data and first presented here. 

d We list combined results for the two components, from Tytler &; Buries (1997) and Buries Sz Tytler 
(1998b). 

e O'Meara et al. 2001. 

'This paper. 

E Pettini & Bowen 2001. 

h Discussed in the the appendix of this paper. 

'First analyzed by D'Odorico et al. 2001. We quote results from Levshakov et al. 2002. 
j From Kirkman et al. 2000. 
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Table 6 

COLUMN DENSITIES AND METAL ABUNDANCES WHERE D/H IS MEASURED 



QSO 


log N H i 


Element 


Abundance 




(cm -2 ) 


Q 


[a/R] 


PKS 1937-1009 a 


17.86 ±0.02 


Si 


-2.7,-1.9 


Q1009+299 b 


17.39 ±0.06 


Si 


-2.4,-2.7 


HS 0105+1619 c 


19.422 ±0.009 


O d 


-1.73 


Q1243+3047 6 


19.73 ±0.04 


o d 


-2.79 ±0.05 


Q2206-199 f 


20.436 ±0.008 


Si 


-2.23 g 


Q0347-3819 h 


20.626 ±0.005 


Si 


-r 


Q0130-403J 


16.66 ±0.02 


Si 


-2.6 



a We list the parameters for each of the two components, where avail- 
able, from Tytler, Fan & Buries (1996); Buries & Tytler (1998a). 

k We list the parameters for each of the two components, where avail- 
able, from Tytler & Buries (1997) and Buries & Tytler (1998b). 

c O'Meara et al. 2001. 

^Using log O/H = —3.31 from Allende Prieto, Lambert & Asphmd 
(2001). 

c This paper. 

f Pettini & Bowen 2001. 
g Prochaska & Wolfe 1997. 
h Levshakov et al. 2002. 
'Discussed in this paper. 
j Kirkman et al. 1999. 



Table 7 

Recent Estimates of the Baryon Density 



Method n b h 2 Ref. 



BBN ■+ 


- D/H 


0.0214 ±0.0020 


this paper 


CMB:- 


BOOMERANG 


021 + 003 


Netterfield et al. 2002 


CMB:- 


DASI 


22+ 004 
u.uzz_ 003 


Pryke et al. 2002 


CMB:- 


MAXIMA-I 


0.033 ±0.013 (95%) 


Stompor et al. 2001 


CMB:- 


CBI 


022+ - 15 


Sievers et al. 2002 


Clusters + SNIa 


n 019+ 007 


Steigman 2002 



33 



f — 1 — 1 — 1 — 1 — I — 1 — 1 — 1 — 1 — I — 1 — 1 — 1 — 1 — I — 1 — 1 — 1 — 1 — I — 1 — 1 — 1 — 1 — I — 1 — 1 — 1 — 1 — I — 1 — h 




4000 5000 6000 7000 8000 9000 

Wavelength (A) 

Fig. 1— The spectra of Q 1243+3047 from the KAST spectrograph (top), HIRES (middle) and ESI (bottom). We show the complete 
wavelength coverage for the Kast and HIRES spectra, but not for the ESI, which extends to 10,000 A. We have applied relative flux calibration 
to all three spectra. The emission lines blend to give a continuously undulating continuum level from 4400 - 5000 A. The vertical marks 
above the Kast and HIRES spectra show the positions of the Lyman series lines in the absorption system at z = 2.526 that gives the D/H 
value. The Lyo absorption line of this system, from which we get the H I column density, is near 4285 A, just to the left of the peak of the 
Lya emission line. We do not plot most pixels, to reduce the file size. 
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Fig. 2. — Expansion of the Kast, ESI and HIRES spectra from Figure 1. The Lya absorption near 4285 A is from the system in which we 
measure D/H. 
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Fig. 3. — Most of the metal absorption lines near z = 2.526. We shifted the Al II 1670 spectrum, which is the only one from ESI, by 
— 16.5 km s _1 to correct a likely error. We see three types of components, grouped by ionization; the low ionization lines represented by O I 
alone, intermediate ions C II, Si II and high ionization C IV and Si IV. 
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Fig. 4. — Voigt profile fits to the lines of the low ionization ions. The line parameters used to generate these profiles are given in Table 3. 
The data is our combined HIRES spectrum, with ltr errors shown just above the zero flux level. Flux is in units of 10~ 16 erg s~ 1 cm -2 A 
— 1 . The O I absorption shows a well defined component 1 at V = 0.0 km s — 1 , and additional absorption in component 2 near +13 km s~ 1 . 
C II and Si II show several components, including a main component near v = km s~ 1 for C II (1335) or 3 km s — 1 for C II (1036) and 
Si II, and a second component near 16 km s _1 . We see no absorption near -82 km s — . 
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Fig. 5. — As Figure 4, but here we show the lines of the high ionization ions. Again, there is no evidence of metal line absorption near -82 
km s _1 . 




Velocity (km sec ) 

Fig. 6. — The blending of the absorption line components D-l, D-2 and H-3 in Lyman series lines 4-9. The two components of the D I 
absorption are centered at —82 km s _1 and —70 km s -1 . H-3 is the deep broad line centered near —40 km s _1 . The data is our summed 
HIRES spectrum. Also shown is the full fit (D-l, D-2, H-3, and all other absorption associated with the system). Note that D-2 is strongly 
blended with both D-l and H-3. 
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The HIRES spectrum of Ly-2 to 8, together with our model of the system, as given in Table 3. 




— The HIRES spectrum 
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of Ly-9 to 18, together with our model of the system, as given in Table 3. 
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Fig. 9. — The b and Ndi values from single component fits to the individual D lines. The only free parameters in the optimizations were 
6(D) and Nrji(D). All other parameters, including v(D)=0, and 6(H-3)=f 6.3km s _1 , were fixed during the fitting process. The displayed 
errors com from the estimated covariance matrix we produced during the optimization process. 
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Fig. 10. — The continuum level used to construct our models of the D absorption. The positions of the Lyman series lines are marked at 
the top of the panel. 




Fig. 11. — The complex continuum shown in Fig. 10 is required to obtain the fit shown here using relatively few Lya forest lines. 
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Fig. 12. — The results of our grid search to measure Ndi- We plot contours of constant x 2 for Ndi against each of the parameters we varied 
in the search. The contours are at the 1,2,3, and 4 sigma levels. 




Fig. 13. — The minimum x 2 value that we found in our grid search, as a function of the total Ndi (cm 2 ) in components 1 and 2 
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Fig. 14. 



Square of the intrinsic width of the lines of ions as a function of 1/mass of the ion in atomic mass units. We plot the b int for 



the main components of the O I, D I and H I (components O-l, D-l and H-l). The solid line connects the O I and H I points and ignores the 
D I. The slope of this line gives a temperature T = 1.1±0.6 X 10 4 K and the intercept gives the turbulent velocity of b tur i, = 5.8±0.6 km s _1 . 
The line predicts b(D — 1) = 11.3 ± 1.8 km s _1 . The observed b(D — 1) = 9.2 ± 0.2 km s is 1.2cr below, and consistent with this prediction. 
The dashed line that is the best fit to O, H and D. The data are also consistent with this fit, which we prefer. Although the b-value of the 
H-l component is not well known, the data shown on this plot provide evidence that D-l is D rather than H. 
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Fig. 15. — Our initial model of the spectrum 4284 A. We show the observed data, our model and the continuum level of our model. The 
jagged horizontal array below 2 units of flux shows five times the absolute value of the difference between our model and the data. The 
squares and circles are the b-spline control points that define the continuum model. The circles are fixed in wavelength and free in flux, while 
the squares are free in both wavelength and flux. The x 2 for this model is 5911, with 4759 degrees of freedom. 
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Fig. 17. — x 2 values as a function of the log Nhi in components 1 and 2. The x 2 values are those returned by the optimization code. The 
results are shown for 3195 restarts. Models containing Lyo forest absorbers with b > 150 km s -1 are excluded from the lower panel. We 
accept log Nhi values that gave x 2 values below the horizontal lines at 6052.3. 
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Fig. 18. — Example of a model with a wide Lya forest absorption line. This model has log Nhi = 19.67 cm -2 and uses the Lya forest 
lines that we show, including the wide, shallow line centered near 4294.5 A with b = 176.7 km s — 1 . The line running through the spectrum is 
the complete model that includes the continuum level, the DLA and the Lya forest lines. The optimizer restarts found a number of models, 
typically with log Nhi < 19.69 cm -2 , which used similar wide Lya forest lines. In these cases, it appears that the optimizer is using wide 
Lya forest absorption to provide opacity that should be provided by the DLA or the continuum. Because such wide Lya forest lines are not 
observed in other studies of the Lya forest, we rejected models with wide Lya forest lines. 
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Fig. 19. — The upper panel shows the normalized in flux profiles of two absorption lines whose column densities differ by log Nhi = 0.04 
cm -2 , the error we quote on the log Njji value for the system that shows D. The upper line has a column density of log Njji = 19.73 cm~ 2 . 
The lower panel displays the difference in flux between these two absorption line profiles. 
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Fig. 20. — As 17, only for the 7898 optimization restarts with three continuum control points in the region between 4260 A and 4325 A. 
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Fig. 21. — Three examples of continuum shapes that were rejected by the filter we applied to the three continuum control point models. 
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Fig. 23. — The continua for the best models we were able to find with log Nhi = f9.69 cm -2 . The continuum that is low at 4290A is from 
the best model with three continuum control points in the region between 4260A and 4325A. The other curve is from the best model with 
two continuum control points. 
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Fig. 26. — Column densities for various metal ions as a function of the hydrogen density nn for a photoionized cloud with log Nhi 
= 19.73 cm -2 , a metal abundance [X/H] = —2.77, and a Jgi2 = 10 — 21 ergs cm -2 s" 1 Hz" 1 sr" 1 . Ionization increases to the right. We show 
solar abundance ratios. The O I column density is insensitive to the ionization. Measured column densities in components 1 & 2 are shown by 
the three large boxes in the upper panel, while four allowed upper limits are shown by arrows. The preferred density is logn^f ~ —1.5 cm" 3 . 
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Fig. 27. — The upper panel shows the neutral fraction of the gas with log Nhi= 19.73 cm -2 as a function of ionization parameter log U 
as returned by Cloudy, assuming a constant Hydrogen density of rin = 0.01 cm -3 . The lower panel shows the size of the absorber along the 
line of sight in kpc, again as a function of log U. In the lower panel, the solid curve that begins at 2 kpc displays the absorber size when 
nu = 0.01 cm~ 3 . The J912 increases to the right along this and other curves with the same shape. We also show curves for nn increased 
(lower on the plot) or decreased by a factor of 10 (dot-dashed lines) or 100 (dashed lines). The steeper curve that begins on the horizontal 
axis at log U = —4.4 shows the size when the J912 = 10 — 21 ergs cm -2 s — 1 Hz -1 sr — 1 . The density increases to the right along this curve. 
The ion column densities indicate log U ~ —2.84. 
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Fig. 28. — Temperature for the gas returned by Cloudy for an absorber with log Nhi= 19.73 cm 2 as a function of the ionization parameter 
log U. 
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Fig. 29. — Cloudy simulations for log Nhi= 16.0cm -2 (left panel), and log Njji= 16.5cm -2 (right) assuming a Hydrogen density nn = 0.001 
cm -2 and a metal abundance of [X/H] = —1.5 as a function of ionization parameter log U. For this figure alone, we enhanced the O and 
Si abundances by 0.3 dex, and lowered N by 0.45 dex to match [N/O] in DLAs and elsewhere (Prochaska et al. 2002; Edmunds, Henry & 
Koppen 2001). The results on the left arc applicable to the component near —40 km s — 1 that we call component 3, while those on the right 
are for the component near +100 km s — 1 . 
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FlG. 30. — Measurements of D/H as a function of the Silicon or Oxygen abundance in the gas. The solid circles are from our group, 
Q2206-199 is from Pettini & Bowen (2001) and the diamond is the local interstellar medium (LISM) measurement (Oliveira et al. 2002). The 
error bars are intended to be la but we suspect that in some cases they have been underestimated. The curves show a closed box model for 
the expected D/H evolution. The solid curve is normalized to the primordial D/H from five QSOs while the dotted curve uses the D/H value 
in the LISM. 
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Fig. 31. — As Fig. 30 but showing D/H as a function of the H I column density. This correlation is unexpected and we believed it is an 
accident. 
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Fig. 32. — Comparison of predicted and measured abundances of four light nuclei as a function of the baryon density. The figure has 
three vertical panels each with a different linear scale. The curves show the abundance ratios predicted for SBBN, from the calculations by 
Buries, Nollett and Turncr(2001). The top curve is the 4 He mass as a fraction of the mass of all baryons, while the three lower curves are the 
number fractions D/H, 3 He/H and 7 Li/H. The vertical widths of the curves show the uncertainties in the predictions. The five boxes show 
measurements, where the vertical extension is the 1<t random error, and the horizontal range is adjusted to overlap the prediction curves. For 
4 Hc the larger box is from Olive, Steigman and Skillman (1997), and the error includes in quadrature the systematic error from Olive and 
Skillman (2001). The smaller 4 He box is from Izotov & Thuan (1998). The D/H box is the mean from five QSOs from this paper. The 3 He 
from Bania, Rood &; Balscr (2002) is an upper limit. The 7 Li is from Ryan et al. (2000). We expect that all the data boxes should overlap 
the vertical band that covers the D/H data. They do not, probably because of systematic errors. 



